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Abstract 
Data from the Automatic Picture Transmissions (APT) of the TIROS-N and NOAA-7 
polar-orbiting satellites has been used to determine cloud-top temperatures and brightness 
for convective storms that produced heavy rain. The study is restricted to the period May 
to September 1979 and 1982 over Alberta and Saskatchewan. 

As a storm passes over a point on the ground, the rain accumulation is a result of 
the combination of the cloud’s speed, size of the rain area, and the rain rate. The results 
presented here suggest that rain accumulation is inversely proportional to temperature, 
due mainly to enhanced rain rates at lower temperatures or a greater temperature 
difference between the cloud top and the tropopause. No correlation was found 
between rainfall and brightness (other than a threshold value of brightness indicative of 
heavy rain), but this was not unexpected as the satellite-measured brightness was not 
corrected for sun or viewing angle. 

Based on the above, a simple rain-estimation technique has been developed. This 
was tested against an independent data set of storms chosen only on the basis of their 
appearance on the infrared satellite imagery, without any knowledge of the associated 
precipitation. The scheme performed well for isolated convective storms, but gave poor 
results for more complicated situations such as merging cells, due mainly to the large 
amount of obscuring cirrus debris. 

The technique presented here can be used as a simple "now-casting” tool for 
isolated storms and should serve well as the basis for a more complicated, and accurate, 


satellite-based, rain-estimation scheme. 
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1. Review and Preview 


1.1 Introduction 


Rainfall is vital to sustain life but too much too quickly can bring floods, destruction 
and death. Precipitation measurements are therefore important for a number of 
endeavors which affect the well-being of mankind. Some of these, discussed by Barrett 
and Martin (1981) in their book "The Use of Satellite Data in Rainfall Monitoring" include: 
rainfall inventories, moisture budgets, water resource monitoring, river monitoring and 
control, flash flood warnings, pest and plague contro! and crop growth programs. 
Unfortunately, the distribution of rainfall is often not well known due to a combination of 
its large fluctuations over short intervals of space and time and a deficient observing 
network. 

There are currently three methods by which rainfall is measured. Most simply, rain 
is collected in a graduated cylinder or bucket, although the design of the bucket may be 
quite sophisticated. These devices (known as rain gauges) are often placed hundreds of 
kilometers apart and inspected only once per day, so that small scale variations of rainfall 
are poorly represented. Near-continuous monitoring of rainfall, in both space and time, Is 
possible with the use of radars capable of detecting precipitation-sized water drops. 
However, only a relatively small area can be covered in this way and, as radars are 
expensive, they are usually located only near population centres. Much effort has 
therefore been directed in recent years to the devising of schemes whereby rain amounts 
can be deduced from satellite imagery of the infared and visible energy from cloud tops. 
Satellites offer the advantage of repeatedly viewing large areas in a relatively short period 
of time (ranging from five minutes to several days). Considerable success has been 
achieved estimating precipitation from convective storms. Synoptic-scale rainfall has 
proved more troublesome, however work in this area continues to progress. An 
excellent review of techniques for inferring rain from satellite images is given in Barrett 


and Martin's 1981 book. Accordingly, only the highlights are presented here. 
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1.2 Review 


1.2.1 Basic Relationships 

The first investigators to relate satellite measured cloud top brightness to 
precipitation were Whitney and Fritz (1961). They were able to show that a small, 
isolated, bright cloud seen on a TIROS-| image was related to precipitation on the ground 
and reasoned that brightness was indicative of heavy convective activity. The major 
conclusion over the next 10 years of study was that the brightest convective clouds were 
the most likely to contain precipitation (Woodley and Sancho, 1971). However, while 
strong radar echoes are associated with bright clouds, the reverse is not always true, and 
the echo area is often much less than the cloud area (Martin and Suomi, 1972). 

Cloud-top temperatures can be evaluated by measuring their infrared radiance and 
assuming black body radiation. Lethbridge (1967), for late spring, found that the 
probability of precipitation within 3 hours after the passage of a cloud increased with 
decreasing cloud-top temperature but was only 50% for temperatures less than -24°C. 
More importantly, this study marked the first time that temperature and brightness 
measurements were combined, which improved results considerably. The 3-hour 
probabilty of rain was only 5% for warm, dark clouds, but rose to 75% for cold, bright 
clouds. It was concluded that brightness was important, as it distinguished between thin 
and thick clouds for low temperatures. A cumulonimbus cell (Cb) could, therefore, be 
identified on a satellite photograph as a cold, bright cloud. 

Cumulonimbus clouds can also be identified by their explosive growth, seen on a 
sequence of satellite images as a rapid increase in the anvil area. Rain-estimation 
techniques which employ this time element are called "life-history" methods by Barrett and 


Martin. 


1.2.2 Life History Methods 

The simplest life-history method is that of Stout et al. (1979). Volumetric rain 
rate is correlated to cloud area and its rate of expansion, cloud area being defined as the 
area enclosed by the 245°K (-28°C) or 200W/m? temperature and brightness contours 


respectively. Coefficients to convert cloud area and its rate of change to rain rates are 
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determined from correlation to radar-derived rain rates. In a test over the southern 
Atlantic, the correlation between radar and satellite estimates of one hour rain amounts 
was 0.84. 

Griffith and Woodley developed a similar technique (described in Griffith et al., 
1978) except that rain intensity is also taken to be a function of the Cb’s stage of 
development. Cloud areas are tracked, classified as being in their growing or shrinking 
stage and correlated to echo area. As radar echoes can be used to infer rain rates, a 
relation can be established between cloud area and rainfall. A different relation exists for 
the growing and shrinking phase, with more rain inferred while the cloud is growing. 
Also, the amount of predicted rain is distributed through the cloud depending on 
temperature, with most of the rain assigned to the coldest areas. The correlation 
coefficients for satellite-estimated to gauge-rain amounts over Florida were as high as 
0.80 for one hour amounts and 0.93 for 6-9 hour amounts. When applied to convective 
storms over the High Plains States (Griffith et al., 1981) rain estimates had to be adjusted 
downward to account for the lower precipitable water content of the air, but otherwise 
the technique worked well. In their 1981 paper, four empirical relationships were 
outlined: 

1. "Rain clouds are those which are as cold or colder than -20°C. 

2 Rainfall is directly proportional to cloud area on any given picture. 

3. Rainfall is inversely proportional to cloud-top temperature. 

4 Rainfall is a function of cloud life cycle, that is, more rain is inferred in the early 
stages of a cloud’s history than in the later stages.” 

Augustine et al. (1981) evaluated the Griffith and Woodley technique over 
southern Florida and found that the correlation of satellite versus gauge rain amounts 
varied drastically for changes in study area and duration. For example, the correlation 
coefficient for an area of 500 km? and one hour was only 0.2, while the maximum 
correlation of 0.9 was obtained for 9350 km? and 6 hours. 

Scofield and Oliver (1977) devised a method by which rainfall estimates could be 
made directly from half-hourly GOES satellite images. Cloud-top temperature is used to 
determine base 30-minute rain rates. Enhanced rain fall is assigned to clouds according to 


features seen on the satellite photographs, such as: rapidly expanding anvils, merging 
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clouds, merging cloud lines and “overshooting tops”, ie. tops of cells that appear highly 
textured on the visible image and which have likely penetrated the tropopause. This 
technique was developed to assist in the estimation of rain amounts associated with 
extreme weather events. It has been succesfully applied to storms which caused flash 
floods over eastern United States (Schneider and Scofield, 1980; Scofield, 1981). 
However, similar to the Griffith and Woodley method, additional factors needed to be 
introduced when the technique was applied to different geographical regions. For 
example, a precipitable water factor was required so that rain amounts over the Plain 
States would not be over-estimated. The opposite problem was encountered over the 
southern states where slow-moving, steady-state storms gave more rain than expected 


(Scofield and Oliver, 1980). 


1.2.3 Bi-Spectral Methods 

Lovejoy and Austin (1979a) introduced a bi-spectral method, that is, one which 
combines infrared and visible-range information. Satellite-measured cloud-top brightness 
and temperature values are assigned to grid points at which the occurrence of rain or 
no-rain is deduced from radar data. Threshold values of temperature and brightness are 
then used to infer rain or no-rain, and varied in order to minimize a "loss function” involving 
the incorrect inferences of: 
ie no rain assigned, rain observed, and 
23 rain assigned, no rain observed. 
In this manner optimum values of temperature and brightness which lead to the best 
estimate of rain areas are determined. 

The next step taken by Lovejoy and Austin (1979b) was to estimate rain amounts. 
They found that rain areas were highly correlated to rain amounts, but weakly correlated to 
rain rates. An average rain rate (as determined by radar) was used to estimate rain 
amounts for convective storms over Montreal and the Global Atmospheric Research 
Programme (GARP) - Atlantic Tropical Experiment (GATE) area. The RMS difference 
between radar and satellite-inferred one-hour amounts were 40% and 20% respectively. 
Using a single image and a single rain rate, it was estimated that the RMS difference for 


one hour convective rain amounts over Montreal would be 49% for satellite versus gauge. 
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It was concluded that the most revealing information concerning rain amount is contained 
in the area of cloud cover, and that a single image provides little information on rain 
intensity. As noted by Barrett and Martin, this runs counter to the claims of the Griffith 


and Woodley or Scofield and Oliver technique. 


1.2.4 Other Work 

Wiley and Laitsch (1983) investigated the effect of varying the amount of satellite 
data on the accuracy of rain estimates. Rain accumulations were estimated for the month 
of August 1979, over the Great Plains, using various frequencies of the number of 
satellite pictures per day (ranging from 2 to 24). Rain rates of 2.5mm /6hr were assigned 
to cloud areas deemed to be raining by a meteorologist, with additional rainfall ascribed to 
storms that appeared more intense. It was found that the correlation of monthly rain 
amounts for satellite versus gauge was 0.51 for 2 pictures/day and 0.59 for 24 
pictures / day. 

Adler and Fenn (1979) used the rate of change of cloud-top temperature (dT / dt), 
measured at 5 minute intervals, to infer updraft velocities. They showed that larger 
updrafts are associated with more severe storms. In 1981, Negri and Adler presented 
examples which showed that the clouds with the coldest tops (Tmin) had both the highest 
radar reflectivity and rainfall rate. They concluded that dT/dt and Tmin were correlated to 
precipitation rate for well-defined Cb’s, as seen on Satellite imagery. Adler et al. (1981) 
presented a physical model that related rain rate to Tmin, along with supporting 
case-studies. Using much of the information contained in the above work, Adler and 
Markus (1983) introduced an index, based on rate of change of cloud area and 
temperature, that could be used to predict severe weather. In hind-casting examples, 
tornadic outbreaks could be predicted as much as 30 minutes in advance. 

Temperature structure, rather than absolute values, may also be used to infer 
severe weather. Recent work (see Hymsfield et al., 1983, for example) has shown that a 
"V" shaped structure is often associated with severe weather, and indicative of a storm 


that has penetrated well into the stratosphere. 
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1.2.5 Physical Relationships 

The rain-estimation techniques described above give only empirical relationships. 
The underlying physical basis is that a cold, bright cloud implies a thick cloud, and that the 
probability and intensity of precipitation increases as thickness increases. That 
precipitation is related to cloud thickness is well known. Rain is usually present in clouds 
thicker than 3,000 meters or those with temperatures below -10°C to -20°C, and rarely 
falls from clouds less than 700 metres thick. It is necessary, then, to relate temperature 
and brightness to thickness. Temperature is considered first. 

Absorption (and re-emittance) of infrared radiation in the 8 to 12 ~“m range is 
strong for liquid water drops but weak for water vapour. The cloud top is, therefore, 
assumed to be a black body so that its radiation is proportional to the fourth power of 
temperature. Measurements from space can then be used to determine the cloud-top 
radiation temperature which is usually equated to the actual cloud-top temperature. From 
knowledge of the atmospheric lapse rate, the height of the cloud top can be calculated. 
For thin clouds, the radiation measured by the satellite's radiometer is a result of some 
layer through the cloud (or possibly even some radiation from the ground) so that the 
radiation temperature is higher than the actual temperature. However, in these cases the 
clouds are usually so thin that precipitation is highly unlikely. More difficult, and of more 
concern here, is the problem of differentiating between a precipitating deep convective 
cell and its associated cirrus anvil, both of which may be radiometrically thick and 
therefore have the same cloud-top temperature. This can sometimes be done using 
brightness in conjuction with temperature. There is little doubt, however, that the 
occurrence of precipitation can be correlated to cloud-top temperatures. This is the 
principal hypothesis on which are built many of the rain-estimation techniques described 
above. 

The interpretation of brightness measurements is much more difficult than that for 
temperature and only a brief over-view is given here. The following discussion is taken 
from Chapter 7 in Barrett and Martin (198 1). 

Studies have shown that brightness is highly correlated to thickness, especially for 
stratus clouds over water. However, clouds become optically thick at approximately 


1km, so the association weakens as depth increases. Brightness from various cloud 
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shapes have been modelled, such as infinite planes, cubes or cuboids, and towers rising 
out of slabs. As well as thickness, the width of the cloud is very important. Reflectance 
from the top of a cube is only one-half that of a plane, but increases as width increases. 
This suggests that the observed increase in brightness of a growing Cb is caused mainly by 
the increase in width. Unexpectedly, turrets act to decrease brightness, as viewed from 
above. The fractional light reflected from the top decreases, but is compensated by an 
increase in radiation from the sides of the turret. 

Generally, brightness decreases as the angle of illumination increases (using O° as 
overhead illumination). However, the rate of decrease depends upon the cloud’s optical 
depth and shape. Also, the brightness, as measured by the satellite's radiometer 
(radiance), depends upon the angle of view. Radiance may increase when the satellite Is in 
the anti-solar direction and the angle of illumination is small. 

At any time, then, the brightness depends on a combination of cloud thickness and 


geometry, angle of illumination and angle of view. 


1.3 Purpose of Present Study 

Rain estimation from satellite data for the Prairie Provinces is in its infancy. The 
only investigation to date was that of Wieler (1981) who, based on comparison with radar 
data, determined threshold values of temperature and brightness (-50°C, 145) that would 
delineate areas of heavy rain. As this study represents only the first small steps, the 
present goals are rather modest. They are: to determine whether precipitation amounts 
or intensities from convective storms are correlated to cloud-top temperature and 
brightness and, if so, is the relationship strong enough to become the basis of a 
rain-estimation technique. 

To this end, storms which gave heavy rainfall (>1O0mm) over Alberta and 
Saskatchewan have been investigated. Cloud-top temperatures and brightness were 
obtained from measurements made by the radiometer on-board the TIROS-N and NOAA-7 
polar-orbiting satellites and recorded at the University of Alberta satellite receiving station. 
Data are recorded for two “passes” per day when the satellite is over Western Canada. 
Passes are separated by about 100 minutes. Rain data consist of amounts and intensities 


supplied by the surface observation network of the Atmospheric Environment Service of 
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Canada (AES). 


Except on a few rare occasions, radar data have not been employed. As this is a 
preliminary investigation, it was thought preferable to keep the analysis of data as simple 
as possible. Also, the study is designed so that case-studies can be quickly and easily 
added at some future date. Radar data do not allow either of the above. Firstly, the 
conversion from radar echoes to rain accumulation maps is very time consuming-almost an 
order of magnitude greater than real time at its worst. Secondly, the rain amounts are 
given only in large ranges (eg 10-25mm) too coarse for use here. Admittedly, this 
difficulty can be overcome, but the data and computer programs are owned by the Alberta 
Research Council, and changes to programs are not always feasible. In fact, at the time 
this report was written, the rain-accumulation program was not completely operational. 
Also, the radar data would have to be re-mapped onto the satellite image (or vice-versa) 
and this is not a trivial task. Finally, and most important, an empirical relationship is used to 
convert radar echoes to rainfall, and the conversion is not necessarily adequate under all 
conditions. The consideration of all possible errors in the radar-derived rain amounts is 
beyond the scope of the present study. 

However, it is undeniable that the inclusion of radar data would be highly desirable, 
and is the next logical step to take. As will be seen, some problems were encountered 
which would not have occurred had radar data been used. Since the number of ground 
observing sites is quite small, the sample size is also small. Moreover, confusion can arise 
if the time of the satellite pass is much different than the time rain was observed. As 
always, when simplifications are made, some compromises are inevitable. 

Chapter 2 presents a description of the satellite and its radiometer. The method 
by which storms were selected and analyzed is given in chapter 3, followed by a 
discussion of each storm in chapter 4. This information is used to construct the 
rain-estimation scheme presented in chapter 5. The application of this scheme to infer 
rain amounts from satellite images is given in chapter 6. Finally, conclusions and 


suggested future work is discussed in chapter 7. 


“% rey ‘oieiia wit tour 
rae Magma RTS sirit ra al J pape areiltaucien nha iene idl 
ois airsions rieidre .Vibhaees Gehow ire an 80) TY 2a RENG eal 
mi guna .ehart sa rot S279. 20) (mor SOF gal e8gors Spral = em 
fi) al arty vel barWwe S16. einagite ler i otgy any furs 20998 YO BCD Yh 
dered wd Toad  .witage? tysets ara ts Binge. of sete") bre JonoowtoY 
Jaeeleege Vietelamoas soe. wey, Mi @e2: OH Cems an? epiiw aaws 
igntew aol 10) Syettd Stee? OFF te Hay B's] SOO! Aven biowy we ment ve “J 
| OT DE BLIP NGI IBA IBeinana os Inara, reatbne yilan slain » Can aut t o 7 
the, “waar. ARON ViMeOW IER fen al OtURAd acf-bas Aeidise at aos Tea aNtID 
mY afeuorm, is DeWoeb~ nbs aft: Aiton attlecoc dn to nolmanignog att  @ = 
. vous Mestre eat 16 sanoh ert ap 
se Veg tigi, a Diy iy MBS WHA? OO Uw i slogan 2h) sevewod 
Da VAIO UT a) ky, arraloga*s ornide ose ig Hy e* algy v gate lsmiboy 2xan ort 2) BAB 
lia ts Satie od aoe .eeu MGad eh ieher Lar he rani BF ad ton hivow risittw 
erie cma MO\AMIAOS BVaRIOM lam OO 2 ORE Bite ee AP Te arp a garia priwsemde. 
2 Sew ede apw thar perp eri! nai ing Fils Ase) eek yiledse writ Toone oN 
aiteilhan 6jo 08rd anhd Hrind oso Ae debscdlane ned ByeNle 
bottan et!  TAeINSiie 20 Gos Mtietng ec? 1 ROU TSS 6 anapeng S vetygtd. | 
& yd wewello? © SBiqedo mi reve 2! basylund Dre “nefsdlee -orew errors nointw ya 
em Jawwenes of beau c WieiiertioiN ait > iIieeto “th miata gee ic noel” 
eho Of Serge ait IO MoMtsa/ciys atT 7 IEtaErID Of RoINUAdS] Setbrog OMEAeRT 7 
bra. eo@thde@nos yee .8 xainde> om need. 4 209K SOheipe aint 2 ' 
\ ite Mi beeedalp zt ew SUtut 


a 


2. The Satellite 
2.1 Orbit and Sensors 


2.1.1 Description 

Cloud-top temperature and brightness was obtained from the Automatic Picture 
Transmissions (APT) data recorded at the University of Alberta from the TIROS-N and 
NOAA-7 polar orbiting satellites. Both satellites are part of the TIROS series and differ 
only in the calibration constants required to convert radiance to temperature (as will be 
discussed shortly). A complete discussion of the satellite and its components is given by 
Schwalb (1978). Consequently, only a brief description is given here. 

The satellite is placed in a sun-synchronous orbit at an altitude of approximately 
865 km. Sun synchronous means that the satellite passes over any given latitude at the 
same solar time each day. The orbital parameters are given in figure 2.1 which also 
displays a typical path of the satellite. The Earth is viewed by an on-board scanning 
radiometer, known as a Advanced Very High Resolution Radiometer (AVHRR), by means of 
a mirror rotating at 360 RPM perpendicular to the satellite's path. Data is obtained across 
the entire horizon but retained only for 55° on either side of the mirror’s nadir. Beyond 
this range, limb darkening and foreshortening limit the usefulness of the data. This 
coverage corresponds to 2997 km on the ground at the equator and extends to 
approximately 55° of latitude at 50° north. As the satellite regresses 25° per orbit, there 
is a large amount of overlap between two successive images. In this manner, Alberta is in 
the satellite's field of view for at least two passes in the afternoon. As the mirror scans 
the Earth data are collected at the rate of 66.5K words per second. However, of this, 
only every third scan line and only two of five frequency bands (one visible and one 
infrared) are fed to the APT transmitter. Furthermore, the data along a scan line are 
averaged to give a near constant resolution of 4km. As only every third scan line is sent, 


the between-line resolution !s: 


(6368 Km)(27)/(102 min)(120 scans min?!) ——Eoucmkt 
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PM ASCENDING NODE ORBIT 


Orbital Period 102 minutes 
Regression 25.5° /orbit west 
Inclination Angle G2 
Height Above Ground 865 km 

Figure 2.1 Orbital parameters for the Tiros-N 


and NOAA-7 satellites. 
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The APT signal is amplitude modulated at an equivalent digital rate of 4160 
words/sec (or 2080 words/line at 2 lines/sec). The information consists of 9309 words 
each of visible and infrared data with the remaining 262 word consisting of 
synchronization pulses and calibration data. This is digitized at the University of Alberta at 
3200 Hz (1600 words/line) so that the 909 words are “squeezed” into only 700 words. 
This reduces the 4km resolution along a scan line to 5.2 km. However, the distance 
between these "new" pixels is only 4.2km. Consequently, some information is smeared 
from one pixel to the next. 

The infrared and visible radiation measured by the radiometer are represented by 
digital counts between 0-255. High counts correspond to high reflectivity and low 
temperature for the visible and infrared bands respectively. Using calibration data from 
the radiometer and the fact that the radiometer output (in digital counts) is very close to 
linear makes it possible to construct a calibration curve which converts counts, as 
recorded by a ground receiving station, to temperature. Brightness counts can be 
converted to reflectance values measured in watts/m’?. However, as outlined in section 
1.2.5, the interpretation of these values is more difficult than for the infrared. As no 
attempt was made to normalize the brightness for angle of illumination or viewing angle, 


the values were left as digital counts. 


2.1.2 Calibration of Infrared Data 

Calibration data are sent from the satellite with each scan line. There are 16 pieces 
of calibration data, each repeated over 8 lines, so that a total of 128 lines are needed to 
receive all the information. This is known as a telemetry frame, as is shown in figure 2.2. 
Of these, only boxes 1-13 and 15 are of concern here. 

The first 9 boxes consist of “wedge” values. Each wedge is assigned a known 
digital count before the satellite is launched. By comparing digital wedge values received 
on the ground to the expected values, it is possible to relate ground-received digital 
counts for cloud temperature and brightness to satellite measured counts. The two can 
differ slightly because of the digital-to-analogue and back to digital conversion of the APT 
signal. Boxes 10-13 are digital counts from the radiometer viewing four platinum 


resistance thermometers (PRTs) embedded in its housing, which is designed to be a black 
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body. Box 15 contains counts for the sensor viewing outer space. To obtain a 

relationship between received digital counts and temperature (for the infrared information, 

of course) the following procedure is employed: 

Le A relationship is obtained between known wedge values and received wedge 
values, using a curve-fitting routine from the Interactive Michigan Statistical 


Library (IMSL) computer package. 


2) Using the curve established above, PRT received counts are converted to satellite 
counts. 
3: Using expressions given by Lauritson et al. (1979) and listed in appendix A, 


satellite PRT counts are converted to temperature. The four temperature 
- estimates are then averaged. 

4. The black-body radiation associated with the above temperature is calculated 
using an approximation of the integral of the Planck function and the spectral 
response of the radiometer (see appendix A). This gives one point for the 
radiance versus digital count calibration curve. 

oe The above point and the counts associated with the radiance of space (box 15) 
are used to construct a straight line graph relating radiance to digital counts. As 
the response of the radiometer is not exactly linear, a small correction is applied 
(see Appendix A). 

(o), Radiance values for each temperature, in 2°C increments from 27°C to -73°C, are 


then calculated using the same routine as in step 4. 


ae Combining (5) and (6), a curve of temperature versus satellite digital counts is 
constructed. 
8. Using curves from (1) and (7), a calibration curve of temperature versus received 


digital counts is made. 
Sh, Finally, received digital counts are converted to temperature, using the curve 
established above and interpolating as necessary between the 2°C-interval 


temperature values. 
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2.1.3 Accuracy of Temperature Calibration 

Two different satellites, both of the TIROS series, are used in this study, each 
having its own calibration curve. These curves, derived as in step (8) above, and averaged 
over 12 passes for NOAA-7, and 18 passes for TIROS-N, are given in Table 2.1. The 
uncertainty quoted is +1 standard deviation and is seen to be of order one 1 count. This 
suggests an uncertainty of +1 count in the digitizing process. Using this, columns 3 and 5 
of Table 2.1 give the uncertainty in temperature, which ranges from 2-3°C at low 
temperatures to 0.5° at high temperatures. It is still necessary, however, to compute a 
separate calibration curve for each orbit. Variations of up to 10°C, for a given count, are 
possible from one pass to the next. 

There are other factors which can contribute to uncertainties in the temperature. 
Even with the non-linearity correction introduced in step (5) above, errors of +1°C at 
300°K (22°C) and -1°K at 200°K (-73°C) can still occur (Lauritson et al., 1979). The sign of 
the error is for actual temperature minus derived temperature. 

Errors can also be introduced if clouds are thin, in which case the radiation 
received by the satellite is from some other layer of the atmosphere, and not from the 
cloud top itself. The derived temperature would be an average temperature through the 
layer. As mentioned in section 1.2.5 this should not be a problem here, as thin clouds will 
likely not be precipitating. A similar condition may exist in a moist, but cloud-free 
atmosphere. Here, some radiation is emitted by the atmosphere as well as the Earth's 
surface. The derived surface temperature will therefore be less than the true 
temperature. As noted by Wieler (1981), Curtis and Rao (1969) estimate that, in the 
Tropics, the magnitude of this effect can be as large as 4°C. As the atmosphere in Alberta 
is drier than that of the Tropics a realistic estimate of the error would be approximately 
2°C. Again, this is not of great importance here, as cloud-free areas are of little interest 
in the present work. 

Combining the various errors, then, the uncertainty in cloud-top temperatures is 


likely to be +3 - 3.5°K. 
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2.2 Gridding of Images 


2.2.1 Resolution of Image 

As mentioned previously, the resolution between and within scan lines is 3.3km by 
5.2km respectively. The smallest area of information, or pixel, represents an area of 
approximately 17 km?. By comparison, resolution at the satellite subpoint before 
processing for APT transmission is 1.1km by 0.8km. The area of a small cumulonibus cell 
may be less than 50 km? or only 3 pixels on the APT image. Clouds of this size are too 
small for the temperature and brightness structure to be properly resolved on the APT 
images. For such cases, where possible, images with full resolution from the High 
Resolution Picture Transmissions (HRPT) of satellite have been obtained from AES. These 
pictures contain a grey scale enhancement curve so that the temperature of any pixel can 
be determined within a certain temperature range. While this does not allow accurate 
interpretation of the temperature, comparison to the APT derived value will give some 


indication if the resolution is deficient. This is discussed in greater detail in chapter 3. 


2.2.2 Gridding Accuracy 

in order to relate satellite measurements to rainfall observed on the ground, it Is 
necessary to locate the position of the rain gauge on the satellite image. This requires 
transforming from a coordinate system described by latitude and longtitude to one 
defined by a scan line number and position along that scan, a procedure which will be 
called "gridding” here. The equations to perform this task are given in great detail by 
Reinelt et al (1975). A summary of the essential equations can be found in Greaves (1981) 
or Wieler (1981) and are repeated here in appendix B. The input required for the 
equations are the orbital parameters given in Figure 2.1 plus the time and longtitude at 
which the satellite last crossed the equator. Basically, the gridding procedure is as 
follows. The time of acquisition of the satellite's signal is recorded and subtracted from 
the time the satellite last crossed the equator. Using the satellite's orbital period and 
regression rate, it is then possible to determine the satellite's position at the time of signal 
acquisition. The same orbital parameters, along with the scan rate (120 scans/min) can 


then be used to calculate the position of any pixel or, as done here, to transform from a 
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latitude-longtitude position to one described by the number of the scan lines after 
acquisition of signal, and the position along that line. 

mele each satellite pass, "grid marks” are applied directly onto the satellite image at 
calculated latitude-longtitude intersections. The positional accuracy of these marks is 
determined by comparing the satellite image to an overlay of prominent landmarks (e.g. 
large lakes) and the known position of latitude-longtitude intersections. The overlay is 
superimposed on the satellite picture and adjusted so that landmarks visible on the satellite 
picture are aligned with those on the overlay. The gridding is accurate if the grid marks 
are co-located with the latitude/longtitude intersections. If not, the gridding is in error 
and can be corrected, based on the magnitude and direction of mis-alignment of the grid 
marks. After correcting for any error, the position of the rain gauge on the satellite image 
can be calculated. When the gridding is poor, the usual reason is uncertainty in the time of 
equator crossing which has been in error as much as 10 seconds (20 scan lines). The 
error normally varies slowly from day-to-day, so rather than applying a large correction 
each day, an approximate correction is applied to the equator crossing time before 
gridding is preformed, and updated as necessary. The amount and need for updating is 
determined from a daily check of the gridding accuracy as described above. It is 
fortunate that timing errors vary only slowly, for when no land marks are visible it is 
impossible to check the accuracy of the grid marks. 

The time of satellite equator crossing is, at best, accurate to only the nearest 


second so that the corresponding uncertainty in distance Is: 


(+0.5 sec)(27) (6368 Km) _ 
(102 min)(60 sec/min) — 3.3 Km 


This is, of course, the distance between scan lines. However, an informal survey of lab 
technicians who perform the gridding operation on a regular basis revealed that the 
smallest discernable shift between grid marks and the overlay latitude-longtitude 
intersections is two scan lines, corresponding to a distance of 1.4mm on the satellite 
images produced at the University of Alberta. A shift of one scan-line is undetectable, 
while a three-line shift is quite noticable. Even so, it is difficult to maintain an accuracy of 


two scan lines across the entire picture. A good match in one area might be made at the 
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expense of a poor match elsewhere. An uncertainty of +3 scan lines (10km) in gridding 
has therefore been adopted. Strictly, this uncertainty should prevail only in a direction 
perpendicular to the scan line. Errors in equator crossing time do not affect the position 
along a scan line. While uncertainties in this direction can be incurred due to errors in the 
longtitude of satellite equator crossing, they are generally quite small and result in an error 
less than the distance between pixels. However, for the sake of simplicity, the 
uncertainty in the position of the rain gauge is assumed to be described by a circle of 
radius 10km. 

Finally, in order to deduce cloud motion, two sequential images must be overlaid 
This can be done using landmarks or the grid marks as a guide. However, due to the 
rotation of the Earth beneath the satellite and the Earth's curvature, the latitude-longtitude 
intersections do not form right angles. This produces a distortion across the image so 
that it is not possible to overlay two images with complete accuracy. If the area of 
interest is small, then the error is usually less than 10km. However, if the area is large, 
uncertainties of 30km are possible. This is of crucial concern when trying to map 


possible rain areas along the path of a convective cell, as will be seen in Chapter 6. 
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3. Selection and Analyses of Data 


3.1 Introduction 

As was stated at the outset, the objective of this study is to correlate 
precipitation, as measured on the ground, with satellite-observed cloud-top temperature 
and brightness of a convective storm. In this chapter the criteria by which case studies 
were selected, sources of data and their limitations, as well as the method used to analyze 


the data are discussed. 
3.2 Selection of Sample 


3.2.1 Required Rain Amounts 

Case-studies were selected on the basis of observed point rainfall, the only 
discriminator being a lower limit on the hourly amount that must be exceeded. A rather 
high limit was chosen for a variety of reasons. A small amount of rain could be due to 
either a weak convective cell passing directly overhead or the result of a more intense 
storm that just grazes the observing site. As large rainfalls are relatively uncommon, 
when they are observed it is more likely that the cell passed directly overhead and the 
observed rain amounts are representative of the storm's maximum output. Also, storms 
which produce large rainfalls are of great general interest as they can, in extreme cases, 
cause local flooding or crop damage. Often, a single storm can contribute over 25% of 
the total monthly rainfall at a given location. Finally, if too low a limit is placed on the rain 
amount, the sample size becomes unreasonably large. AES, for the purposes of weather 
reporting, defines heavy rain as arain rate of > 7.5mm/hr. This, of course, refers to rain 
intensity. There are no guidelines as to what accumulation constitues a “large” amount so, 
somewhat arbitrarily, 10mm has been adopted here. There is no objective basis for this 
choice. However, amounts of 5mm or less are fairly common and may pass without 
notice, while 20mm within an hour is a rare event and generally receives considerable 
attention. 

The requirement for a storm to be included in the sample, then, was a report of > 


10mm of rain in one hour or less. A further requirement was that the rain must have 
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occurred within one hour of the satellite pass. 


3.2.2 Area of Study 


It is desirable to confine the study area to regions which fall under the same 
weather regime and climatic controls. Precipitation amounts may be dependent on 
geographical location without being reflected in the cloud-top temperature or brightness. 
An example is the lower rainfall over the Plain States as compared to the eastern United 
States due to the lower water vapour content of the air. Relatively uncomplicated 
topography is also required. Mountains and valleys contribute to mechanisms which can 
result in great variations of local weather and precipitation patterns over small distances. 

The above conditions refer to the need for the study area to be homogeneous. As 
the area becomes large, there is a greater risk that some unknown variable will affect 
precipitation amounts. Homogeneity is of great importance, but, from a practical 
viewpoint, if the area is made too small then the sample becomes too small to be of any 
value. This is especially true here, as there are only a few observation sites distributed 
over a very large area. Consequently, the study area is of considerable size, namely 
Alberta and Saskatchewan. These two provinces are characterized by flatland or rolling 
hills in the south, a mixture of forests and lakes in the north and the Rocky Mountains to 
the west. However, even within a relatively homogeneous region the climate may vary 
considerably. Southern Saskatchewan, for example, is much drier than most of Alberta, 
and the area near Red Deer, in central Alberta, is well known for its frequency of 
hail-storms. 

British Columbia was excluded because of its complex terrain. Manitoba was not 
included as it can come under the influence of warm, moist air from the south, a condition 


which is uncommon in Saskatchewan and very rare for Alberta (Hare and Thomas 1974). 
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3.3 Sources of Rain Data 


3.3.1 Tipping Bucket Rain Gauge 

Ideally, a running total of rain accumulations over a very short time period (i.e. of 
order one minute) is best suited for comparison to satellite pictures. This would allow 
calculations of rain intensities while the satellite was overhead, as well as accumulations 
over the time period between passes. At selected observing sites this information is 
recorded using a device known as a tipping bucket rain gauge (TBRG). The gauge is 
connected to a rotating drum (one revolution in 24 hours) on which a tick is placed for 
each 0.2mm of rain that falls into the bucket. Unfortunately, these records are not easily 
accessible. Instead, hourly amounts and maximum accumulations for time periods ranging 
from 5 minutes to 12 hours are calculated from the TBRG chart, and monthly summaries 
are made available to the public on a form known as TAB31, an example of which is given 
in figure 3.1. The right hand side of the chart gives the hourly accumulations in tenths of 
millimeters. The left hand side gives maximum accumulation over various time periods 
which can then be converted to rain intensities in mm/hr. However, at best, the time that 
these values occurred is known to only the nearest one-half hour, and then only if there 
was no more than one period of heavy rain that day and it lasted less than one hour. This 
would be the case, for example, if just one convective storm passed over the rain gauge 
on any given day. Fortunately, this is a relatively common occurrence so that the left hand 
side of the TAB31 chart can be used to estimate both amount and duration of 
precipitation. 

Column 3 of the TAB31 chart is the ratio of the 24-hour rain accumulation as 
measured by a standard gauge and the TBRG. This is used to correct the TBRG amounts 
so that there is only one estimate of 24-hour rainfall. When the factor is outside the 
range 0.5 - 1.5, the amounts are left uncorrected and flagged as such. Normally, this 
occurs only when rainfall totals are very light. 

There is no guarantee, though, that the rain accumulation has been accurately 
measured by the standard gauge. Under windy conditions, when rain is blown away from 
the vertical, the gauge amount can be a gross under-estimate of the true rainfall. Even 


when it is not windy, large errors may be incurred. Barrett and Martin (1981) give results 
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Figure 3.1 
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Monthly Report (TAB3 1) form. 
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from various investigators which indicates that the difference of measured rain amounts 
from co-located gauges is in the range 5 - 20% but can be as large as 30%. When hail 
accompanies the rain, errors are probably larger. Hail cannot enter the TBRG but can 
easily fall into the standard gauge. The following procedure is used to determine the 
amount of liquid precipitation. Hail is collected from an undisturbed area on the ground 
(equal in size to the opening of the gauge) and its water equivalent measured. This is then 
subtracted from the total amount of precipitation (hail plus rain) in the standard gauge. 
Under these conditions, uncertainties in rain amounts of 20-30% are expected. 

While such large uncertainties are disquieting, they do not present a major 
drawback. Rather, they represent a limit on the accuracy of any satellite-inferred rain 
amounts. The ultimate goal to strive for then is for predicted and measured rain amounts 
to agree within 20%. 

The geographical distribution of stations equipped with TBRGs is given in figure 
3.2. Not all stations are operational at the same time as they open and close as economic, 
industrial and social needs dictate. The total of stations given in figure 3.2 is 63, with 31 
in Alberta and 32 in Saskatchewan. Normally, approximately 59 of these are operational 
at any given time. The majority of stations are south of latitude 55 with only two stations 
in the mountains. As aresult, most case-studies are from the southern half of Alberta and 


Saskatchewan. 


3.3.2 Climate Stations 

Daily rain accumulations are available from the AES climate station network shown 
in figure 3.2. The 24 hour period is normally 08:00 to 08:00 local standard time. These 
can be used as an aid in verifying satellite-inferred rain amounts in that they represent an 


upper limit to the predicted amounts. This will be discussed in more detail in chapter 6. 


3.3.3 Surface Weather Reports 

There are a number of observing stations (see figure 3.2) which report weather 
conditions every hour and more often if there is a “significant” change. The significant 
changes of concern here are: begining or end of precipitation, beginning or end of a 


thunderstorm and certain changes of visibility or cloud ceilings. Changes in rain intensity 
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Figure 3.2b Precipitation network over Saskatchewan 
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#- stations equipped with tipping bucket rain gauge 
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need not be reported, but they are often accompanied by significant changes in ceilings or 
visibility. From these reports it is possible to obtain the duration of precipitation and, 
more importantly, the time it occurred accurate to within a few minutes. This is much 
better than the one-half hour provided by the TAB31 chart. Important information can 
also be obtained from the remarks section of the reports, such as movement of 
thunderstorms or location of heavy showers. Finally, cloud amounts may reveal whether 
or not a storm passed directly overhead. For example, a report that showed only a few 


tenths of the sky covered by cloud would indicate that the cell only grazed the station. 


3.4 Selection Effects 

Attention is now turned to some implications induced by the method of selecting 
the sample. 

A typical thunderstorm cell goes through three stages of development: growing, 
mature and dissipating. Most of the rainfall occurs in the mature stage, and very little is 
associated with the dissipating stage. Often, three cells are present in a single convective 
storm, each in a different stage. The lifetime of a single cell may be less than one hour, 
but because new cells are being generated as old ones dissipate, the storm may endure 
for several hours. 

By choosing the sample based on rain amounts > 10mm, the sample assuredly 
consists mainly of mature storms. However, when a rain-estimation technique is applied 
to randomly selected satellite images, storms in all stages of development might be 
represented. A correct prediction of rain for a storm in its mature stage would be a vast 
overestimate should the cell be in the dissipating stage. For any rain-estimation scheme to 
have a chance to succeed it must be able to distinguish between the different stages of 
development. 

Normally, the stage of development can be determined from a sequence of 
satellite images. Scofield and Oliver (1977), for example, characterize a dissipating storm 
by the cloud area shrinking or stationary with time. Unfortunately, with a time resolution 
of 100 minutes, this procedure is not feasible here. 

By necessity, a different approach is used. A growing storm will not be as “cold” 


as a mature storm as the cloud top Is still rising upwards. Similarly, due to the large thin 
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anvil, a dying storm will be less bright than a mature storm. It may be, then, that some 
broad limits of temperature and brightness can be used to delineate a cell in the mature 
stage. This is an imperfect solution at best, but without more frequent satellite pictures, it 
is the only alternative. 

Another difficulty imposed by the selection method is that to infer rain amounts 
successfully requires not only that mature storms are cold, bright clouds but also that all 
cold, bright clouds produce rain. The sample selection procedure excludes the possibilty 
of finding cold bright non-raining cells. To detect these would require selecting cases 
strictly from the satellite images. However, the useable cases would still be restricted to 
those in which the storm passed directly over a tipping bucket rain gauge. Considering the 
small number of gauges, a large number of storms would have to be considered just to 
find these few fortunate events. As a compromise, a small number of days have been 
chosen using only the satellite imagery (the independent sample) to test against any results 


found from cases selected using rain data (the dependent sample). 


3.5 Analyses of Satellite Data 

Convective storms suitable for study are determined from the hourly amounts 
given in the TAB31 charts. For these storms, contour maps of satellite viewed 
temperature and brightness are produced. The next step is to determine whether or not 
the convective core passed directly over the station. Moreover, when the time of the 
satellite image and observation of rain are considerably different, it is necessary to adjust 
the position of the cloud, in order to estimate the temperature and brightness over the 


station at the time rain was observed. 


3.5.1 Position of Raingauge Under Cell 

The relative position of the rain gauge beneath the convective cloud is of great 
importance in determining the expected amount of rain. The most copious rain falls 
beneath or very close to the deep convective core, with lesser amounts associated with 
the anvil. On a satellite image it is often difficult to separate the core from the cirrus 
debris in which it is embedded. There are some clues which make this task easier, as 


outlined by Scofield and Oliver (1977). In most cases the core is usually found near the 
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Figure 3.3 Relation of a typical thunderstorm’s rainfall pattern 
to the cloud-top temperature contours. 
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tightest temperature gradient on the upwind side of the cloud. This is shown 
schematically in figure 3.4. However, if the cell is small or growing so that the anvil has 
not yet spread downwind, then the core is in the middle of the cloud. Another difficulty is 
determining the size of the rain area. When inferring rain from satellite images over 
Montana, Griffith et al. (1981) found that apportioning most of the rain to only 1 - 10% of 
the cloud area best fit the radar-derived rain areas. In the same region, Wiley and Laitsch 
(1983) assigned rain to 25 - 50% of the cloud area. These studies were estimating rain 
amounts, but the situation here is somewhat different. All that is required is some means 
of dividing the cloud area in order to separate the most likely area which contains the 
convective core from the rest of the cloud. Using the two studies noted above as 
guidance, this has been done by dividing the convective cloud into three zones. Zone 1 
comprises 20% of the cloud’s area nearest the upwind temperature gradient, or the 
central area if there is no large elongated anvil. This is presumed to contain the convective 
core, though it !s possible for the core to be much smaller than this area. Zone 2 
represents a transition between the core and the anvil under which rain is expected, even 
if the storm only grazed the station, and consists of 30% of the cloud area. Zone 3 
encompasses the remaining 50% which represents the thin, non-precipitating anvil. 

Case studies are classified according to which cloud zone passed over the station 
in order to remove the dependence of rain amounts on distance from the convective core. 
When two or more cells have merged the brightness can often be used to separate the 


storms before attempting to construct the zones. 


3.5.2 Adjustment of Cloud Position 

When the time of the satellite pass and observation of rain are nearly coincident, it 
is immediately apparent which cloud zone passed over the rain gauge. However, when 
they differ, the motion of the storm needs to be known. Given two sequential satellite 
images this is a simple task, assuming straight-line motion, a reasonable assumption over 
the short intervals considered here. However, when cells are merging or evolving rapidly, 
then it may not be possible to identify which zone passed over the station. For these 
cases, two zones may be quoted with the most likely given first. When only one image !s 


available, storms are assumed to move with the 7O0Omb flow. As this is only an 
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approximation, two possible zones will be given for the majority of these cases. 

The adjustment of the cloud position is also necessary for determining the 
temperature and brightness over the rain gauge at the time rain was observed. If the time 
of rainfall is not accurately known, then it is assumed that the rain is associated with the 
lowest temperatures that passed over the gauge, as determined from the track of the 
storm. In all cases, after adjusting for cloud position, the temperature associated with the 
rain is assumed to be the lowest within 10km of the gauge, to account for position 
uncertainty. This is especially important near areas of tight temperature gradients. The 
positional accuracy here is insufficient to determine if the rain is associated with the 


gradient or merely an artifact of the uncertainty in the gauge’s position. 


3.5.3 Evolution of Cloud 

While the above procedures can be used to correct for the position of the cloud, 
it Cannot account for any evolution that may have taken place in the interval between 
observation of rain and the satellite pass. 

Adler and Fenn (1979) showed that, for severe storms, the cloud-top 
temperatures can decrease faster than 1°K/min. through the range 240°K to 235°K (-33°C 
to -38°C). However, the associated severe weather occurs 15 to 30 minutes after the 
period of maximum growth, when the minimum temperature is reached. From their 
published curves, below approximately -45°C to -50°C temperatures decrease at a much 
slower rate (as low as 0.1°C/min). Accordingly, reasonable temperature estimates (+5°C) 
are possible if the satellite pass is within 30 minutes of rain observation. Furthermore, 
both the Scofield and Oliver, and Griffith and Woodley techniques use 30-minute GOES 
data to infer rain successfully. This would also indicate that, except for a rapidly growing 
cell, the temperature at a given instant is representative of conditions over one-half hour. 
Unfortunately, no parallel study exists for brightness, and its rate of change is unknown. 

Generally, the longer the time between satellite pass and observation of rain, the 
larger the uncertainty of the quoted temperature and brightness values. However, as the 
rate of evolution for any particular storm studied is unknown, the maximum time 


difference for which results become unreliable is also unknown. 
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3.5.4 Additional Sources of Satellite Data 

As was noted earlier, the APT images consist of only every third scan from the 
AVHRR and that the data along each scan are averaged. HRPT images employing all scan 
lines with full resolutions are available from AES. 

A "grey scale” is applied to the HRPT images so that each pixel is represented by a 
certain shade of grey, dependent on its temperature or brightness. The calibration of the 
temperature grey scale is made using a curve similar to that given in Table 2.1. An 
enhancement curve can be constructed so that certain temperature ranges are highlighted, 
but for most of the infrared pictures and all of of the brightness images used here, the 
0-255 range of digital values is divided evenly into 17 shades of grey. This leads to 
increasingly larger ranges of temperature at high digital counts. For example, the -40°C to 
-56°C range is represented by only one shade of grey. 

While such large ranges might appear to make the HRPT images useless, this is far 
from true, due to their superior resolution. The resolution at the satellite subpoint is 
approximately 1.2km? compared to 17km? for APT pictures. Small cloud areas that may 
be poorly resolved on the APT images, especially those with tight temperature gradients, 
may be better represented on the HRPT photographs. If the HRPT and APT temperature 
estimates differ by more than one shade of grey, the APT estimate is disregarded because 


of its lower resolution. As will be seen, this occurred on two occasions. 


3.5.5 Brightness 

As mentioned earlier, brightness values are left as digital counts, and have not been 
corrected for sun angle. As the specific cases examined here are all confined to a 
one-month period only (though not of the same year) at approximately the same time of 
day, the sun angle is approximately constant. However, because brightness values have 


not been “normalized”, comparisons to most other studies cannot be made. 
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3.6 Application of Sample Selection and Analyses 

The previous sections have explained the procedure used to obtain and analyze the 
dependent sample. Given below are the steps taken in the application of these 
procedures. 

The time period for selection of the sample was restricted to May through 
September, as convective activity is rare outside this period. Satellite data were availabe 
for June -Sept. 1979 and May-Sept. 1982. Locations for which > 10mm of precipitation 
in less than one hour were reported were found using the TAB31 charts (charts for 
Saskatchewan had to be ordered specially and those for May 1979, 1982 and June 1982 
were not received). Only those storms for which the precipitation was within one hour of 
the satellite pass were considered. Next, satellite pictures were then inspected to 
determine if the precipitation was the the result of convective activity. This was generally 
not a problem; isolated cells or obvious cells embedded in other clouds were evident in all 
but two cases. 

Temperature and brightness contour maps were then produced and the location of 
the rain gauges placed on the maps. Finally, the cloud position was adjusted as explained 
in section 3.4.2 and the temperature and brightness over the gauge estimated. Where 
other convective clouds, not selected by the procedure given above, were found to lie 
near a rain gauge, the temperature, brightness and precipitation amount were also 
recorded. These additional storms serve two purposes. They may aid in determining the 
temperature and brightness cut-off for which < 10mm is expected or reveal the existence 
of cold, bright, non-precipitating storms. 

Each case is discussed in some detail in the next chapter. Even cases which 
eventually proved not to be useable are included, as they illustrate some of the difficulties 


encountered. 
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4. Case-Studies 


4.1 introduction 

For each case, temperature and brightness contour maps are given. Unless 
otherwise indicated, the scale of all maps is approximately 1 cm = 45 Km. An arrow is 
drawn in the upper-right hand corner of each map pointing northwards. Generally, only 
the O°C and 80 digital count contour, for temperature and brightness respectively, are 
labelled. Shading is used to indicate areas enclosed by temperatures below and brightness 
above some "critical’ value. The value changes from one case to the next and is given in 
each figure. The contour intervals are given in Table 4.1. When available, weather 
observations were obtained and these are listed in appendix C. At the end of the chapter, 


a summary table of all pertinent information is given. 
4.2 Case Studies 


4.2.1 July 07 1979 

Edson received 11.1mm from 20-21 GMT, with the heaviest precipitation 
occurring from 20:05 GMT to 20:40 GMT, as seen from the weather observations given 
in Table C.1. Edson was scanned at 20:30 GMT on the first of two satellite passes over 
Alberta on this afternoon. The temperature and brightness maps are given in Figure 4.1. 

A large storm is seen to lie over western Alberta. As Edson was scanned while 
the heavy rain fell, the temperature and brightness associated with the rain are those 
values directly over the station, namely -53°C and 170 respectively. Lower temperatures 
of -57°C lie just 20km to the northwest. Even though the cloud extends far south of 
Edson, it is very narrow and the area encompassed Is quite small. Consequently, Edson is 
determined to have been passed by Zone 1 of the cloud. 

The second satellite image is at 22:10 GMT. Temperature and brightness maps are 
given in Figures 4.2. It is seen that the storm decayed while moving northeast at 
approximately 50km/hr. Minimum cloud-top temperatures decreased 5°C. No heavy rain 


was being reported at this time. 
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4.2.2 July 10 1979 

The temperature and brightness maps for the 21:40 GMT satellite pass are given in 
Figure 4.3 and 4.4 A number of stations recorded precipitation near this time with both 
Edson and Saskatoon receiving at least 10.0mm. Only one image is available for this day. 
The upper air flow (not shown) was predominantly from the southwest, except over 
central and northern Alberta where it was light and variable. 

Figure 4.4 shows a cloud band with embedded ceils running through central 
Alberta before curving southward over Saskatchewan. A number of isolated cells are 
also evident over Saskatchewan as seen in Figure 4.3 Each area of interest is discussed 
separately below. 

In Alberta, Edson is under Zone 1 at the southern end of a bright cell. The 
temperature and brightness over the raingauge are -49°C and 156, respectively. 
Brightness values in excess of 170 are only 10-20km to the west. 

Edson received 11.3mm of rain from 21-22 GMT. The weather observations, 
given in Table C.2, show that only light rainshowers were present at 21:00 GMT, but 
moderate showers were occurring at 22:00 GMT. While no heavy rain was reported, it is 
assumed to have occurred closer to 22:00 GMT. This is consistent with Edson under 
Zone 1 at the time of the satellite pass. Accordingly, the temperature and brightness over 
Edson from the contour maps have been adopted as being associated with the heavy 
rainfall. It is also noted, from the 2200 GMT report, that the storm was moving to the 
northeast. 

Another large cell lies over central Alberta. At the time of the satellite pass, 
Calgary was just south of this cell with Penhold and Rocky Mountain House east of the 
cloud. Assuming a northeasterly motion, both Penhold and Rocky Mountain House would 
have been passed by Zone 3 of this cloud. It is possible that the storm missed Calgary 
completely. 

As expected, very little rain was observed at these three stations. Calgary did not 
receive any rain from 21-22 GMT, while Rocky Mountain House and Penhold had only 
0.8mm and 2.6mm, respectively from 22-23 GMT. 

In southern Saskatchewan, a complex of three cells lie near Swift Current. The 


high resolution AVHRR visual image shows that there are distinctive breaks between the 
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cells. This is also detectable on the APT images, but is not as obvious. There are two 
TBRGs in Swift Current, one at the airport and another operated by the Canadian Dept. of 
Agriculture. Both lie very close to the breaks between the cells, with the temperature and 
brightness overhead equal to -49°C and 143 for Swift Current A, and -49°C and 153 for 
Swift Current CDA. 

Swift Current CDA received 4.4mm from 20-21 GMT and 3.9mm fell at the 
airport during the following hour. A northeasterly cloud motion would result in the storms 
just grazing the gauges, so that they would have been passed by Zone 2 or 3 of the cloud, 
although Zone 2 seems more likely. Considering the proximity of the cloud over the 
airport, the temperature and brightness values overhead are accepted as representative of 
the conditions prevalent when the rain was observed. As the rain at the agricultural station 
occurred at least 40 minutes before the satellite pass, it will not be considered here. 

The final cloud of interest is the large cell lying over central Saskatchewan. 
Outlook, at the southern end, would have been crossed by Zone 2 or 1, while Saskatoon 
is under Zone 1 and North Battleford is beneath the anvil. The temperature and brightness 
over Saskatoon are -52°C and 165, respectively. The speed and movement of the cell are 
unknown; however, assuming that the storm moved northeast, then the temperature and 
brightness values that would have passed over Outlook are -45°C +5 and approximately 
160. 

Saskatoon received 10mm in just 15 minutes. The weather observations, given in 
Table C.2 show that heavy showers occurred from 21:03 to 21:26 GMT, which is close 
to the time of the satellite image. As the gradients are quite flat, the temperature and 
brightness values over Saskatoon are representative of those associated with the heavy 
showers. 

Outlook received 6.1mm from 20-21 GMT with 5.7mm in only 5 minutes. North 
Battleford recorded 3.7mm from 19-20 GMT. Weather reports for North Battleford, 
given in Table C.2, show that heavy rain was observed, but 50 minutes before the satellite 


pass which is to long to be of use here. 
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As in Figure 4.3 except over Alberta. 


Figure 4.4 
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4.2.3 July 11 1979 

This is one of two cases where the resolution of the APT data was insufficient to 
measure the cloud-top temperature accurately. (see also July 15, 1979). The rain-gauge 
data from around Edmonton confirm that the disturbance treated here is of a very small 
scale. Edmonton Municipal Airport received 39.5mm of rain between 19-22 GMT. 
During the same time period Edmonton Namao, Ellerslie and Edmonton International Airport 
received 8.0mm, 2.1mm and 0.6mm, respectively. From 20-21 GMT, the Municipal 
Airport recorded 17.7mm, while the other three stations received barely more than a 
trace of precipitation. In the next hour (21-22 GMT) a further 10.3mm fell at the 
Municipal Airport. The weather observations, given in Table C.3, show that the Municipal 
Airport reported heavy thundershowers from 21:26 GMT to 21:32 GMT. The movement 
of the shower was toward the east or southeast. Namao reported heavy showers to the 
south at 21:00 GMT and thunder at 21:37 GMT but no precipitation until 22:00 GMT. 
During this time the International Airport reported heavy showers to the north. At 
Ellerslie, the University of Alberta climatological station, only 1.3mm fell from 21-22 GMT 

Attention is now turned to the temperature and brightness maps for the 20:35 
GMT satellite pass, shown in Figure 4.5 and 4.6. No large convective clouds are seen 
over Edmonton although considerable activity is evident to the north. Both Namao and the 
Municipal Airport are under an area of relatively high temperatures, whereas it is known 
that heavy rain was occurring at this time. The nearest clouds to the Municipal Airport lie 
less than 10km to the west, with a cloud top temperature of -27°C. However, close 
inspection of the HRPT infrared image for this day (not shown) reveals a small bright spot 
in the centre of the cloud which is definitely a whiter shade of grey than the surrounding 
cloud. Comparison to the grey scale gives a temperature in the range -40°C to -56°C. 

lt is assumed that this is the storm which dumped heavy rain on the Municipal 
Airport and that positioning errors account for its apparent location west of, rather than 
over Edmonton on the satellite image. Comparing the APT and HRPT data, the derived 
temperature differs by > 13°C, indicating that the cloud is below the limit of resolution of 
the APT image. 

While this case is left in the sample (at least for the time being) it is of limited use 


as the temperature is known only to within +8°C. The maximum brightness of the cloud, 
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from Figure 4.5b, is 165. A visible range HRPT image was not available for this day, and 
so due to lack of any other information the above brightness value is accepted as found. 
Considering the large amount of rain, it would seem safe to assume that the Municipal 


Airport was crossed by Zone 1 of the storm. 
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4.2.4 July 14 1979 

On this day, Kindersley received 8.6mm of rain from 21-22 GMT. Temperature 
and brightness contour maps for the satellite passes at 21:50 GMT and 22:30 GMT are 
given in Figures 4.7 and 4.8, respectively. 

Based on the surface weather reports, given in Table C.4, most of the rain 
occurred near 21:30 GMT. This is 40 minutes after the first pass and 60 minutes before 
the second. Due to this large time lag, it is difficult to discern which cloud passed over 
Kindersley and produced the precipitation. By overlaying the two maps, a general 
southerly motion of clouds in western Saskatchewan is evident, but it is very difficult to 
follow any particular cloud area. The direction of flow at 500mb and at the surface (not 
shown) also supports this movement. 

The disorganized nature of both the temperature and brightness fields belie the 
convective nature of the precipitation (as reported in the surface observations). Even 
though it is known that Cb’s are present, this is far from obvious from Figures 4.7 or 4.8. 
It is possible that they are hidden by a layer of higher cloud. However, if Cb’s are present, 
then their tops must be quite warm, as no particulary cold cells are evident near Kindersley 
on either image; temperatures being in the -30°C to -35°C range. 

The lack of any well defined cells and the large time lag between the satellite 
images and the observation of rain over Kindersley make it impossible to identify the 
rain-producing cloud. Consequently, this case it not used. However, it is important to 
note that cloud tops are quite warm. It is unknown if insufficient APT resolution played a 
role, as HRPT images were not available. The surface observations show that there were 
many brief showers, suggesting that the weather was dominated by many small cells. It is 
possible, but unsubstantiated here, that either some small cell beneath the resolution of the 


APT image or a cell with a short lifetime (or both) caused the heavy rain at Kindersley. 
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As in Figure 4.7 except at 2230 GMT. 


Figure 4.8 
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4.2.5 July 15 1979 

The satellite images at 20:45 GMT and 22:25 GMT show only scattered small cells, 
as seen on temperature and brightness maps given in Figures 4.9 and 4.10. However, 
when the satellite images are superimposed it is evident that one of the small clouds did 
indeed pass over the rain gauge. Outlook received 11.7mm in only 15 minutes, sometime 
between 22 GMT to 23 GMT. 

From Figure 4.11, the minimum cloud top temperature is only -35°C, quite warm 
for a suspected Cb, and the small cloud size raises suspicions regarding possible 
problems with resolution. From the HRPT infrared image for this day (not shown), 
comparison with the grey shade gives a cloud-top temperature in the range -40°C to 
-56°C. In fact, a few co-workers asked to determine the category for this cloud placed it 
in the next coldest range. 

The temperature adopted !s -48°+8°C with a brightness, from fig 4.10 of 1555. 
However, similar to July 11,1979 while included in the sample, the large temperature 


uncertainty limits the usefulness of this case. 
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Figure 4.11 Expanded view of temperature contours near Outlook 
on July 15, 1979 at 2225 GMT. 
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4.2.6 July 29 1979 

The temperature and brightness maps for the 21:40 GMT pass are given in Figure 
4.12 and show Calgary lying under Zone 1 of a large cloud. The temperature and 
brightness values over Calgary are -43°C and 168, respectively. However, the lowest 
temperatures of -48°C are only 10km to the southwest. 

The second satellite pass at 23:00 GMT viewed only the western half of Alberta 
(not shown). A cloud mass with a temperature of -50°C is seen at the eastern edge of the 
image. If, as seems probable, this is the same cloud that was over Calgary, then cloud 
motion is easterly at approximately 45km/hr. 

Calgary received 10.0mm of rain from 21-22 GMT with 8.1mm in 30 minutes. 
Weather reports for this day were not available, but there is little doubt that the cell over 
Calgary on Figure 4.12 produced the rain. In accordance with the discussion given in 
section 3.4.2, the lowest temperature within 10km of the gauge has been adopted as 
associated with the rainfall, namely -48°C. The brightness, however, remains unchanged 


at 168. 
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4.2.7 August 05 1979 

The temperature and brightness maps for the 20:25 GMT pass are given in Figure 
4.13. A large, organized band of cloud with embedded Cb's, lying along a TROWAL, is 
seen to curve northwest-southeast through Alberta and Saskatchewan. Other isolated 
cells cover western Alberta and southern Saskatchewan. 

Cold Lake lies underneath a cold, bright area within the cloud band. The 
temperature and brightness overhead are -57°C and 153, respectively. Temperatures of 
-60°C and brightness values over 160 are only 20km to the west-southwest. 

Cold Lake received 24.0mm from 20-21 GMT, all within 30 minutes. Weather 
reports for this day were not available, so a more accurate timing of the rainfall is not 
possible. For lack of any better information then, the values of temperature and 
brightness over Cold Lake have been adopted as associated with the rain. 

From the temperature and brightness maps for the second pass, at 22:05 GMT, 
(see Figure 4.14) it is seen that the edge of the cloud band has moved approximately 80km 
eastward. The movement of the individual cells, however, appears to be to the 
northwest, but this is not certain. The cloud mass that was over Cold Lake has dissipated 
and merged with the background cloud. Consequently, its direction of motion is unknown 
and Cold Lake may have been passed over by Zone 1 or 2 of the storm. 

Two cells are seen over western Alberta in fiqure 4.15. Edson lies directy under 
Zone 1 of one of these, with a cloud-top temperature of -50°C and brightness of 156. 
From 21-22 GMT, Edson received 14.1mm of precipitation. 

Rocky Mountain House is under Zone 1 of the other cell, with a cloud-top 
temperature overhead of -35°C. The minimum cloud temperature is -40°C. The 
brightness over the station is 155 with values of 160 only 10km to the north. Rocky 
Mountain House recorded just 0.7mm from 21-22 GMT. 

Based on the satellite maps, a cold, very bright storm passed directly over 
Outlook. In fact, the brightness value of 190, seen on Figure 4.14, was the highest value 
encountered for the entire study. Unfortunately, hourly rain amounts are missing for this 
day. For interest’s sake only, the 24 hour accumulation was 14.1mm, though how much 


of this is due to this storm is unknown. 
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4.2.8 August 21 1979 

The first satellite pass was at 21:00 GMT. From the temperature and brightness 
maps, given in Figure 4.16, a small Cb is evident 55km northwest of Lethbridge with a 
cloud-top minimum temperature of -54°C and maximum brightness of 161. 

By the time of the second pass, at 22:40 GMT, considerable changes have 
occurred. The temperature and brightness maps in Figure 4.17 show a very cold cloud 
mass of low brightness covering southern Alberta. The coldest area lies close to the 
position of the Cb noted from the first pass, with a temperature of -64°C and brightness 
of 140. If this is the same cell as on the first image, then it has moved very little indeed. 
While the 500mb flow, as shown in Figure 4.18 is quite slack, it is unusual for a Prairie 
convective storm to remain stationary for over 100 minutes. More likely, the Cb noted 
from the first pass drifted southeast while a new cell developed in its place. The 
temperature and brightness over Lethbridge at 22:40 GMT are -55°C and 127, 
respectively while the highest brightness values in southern Alberta (155) are 15km south 
of Lethbridge. It is suggested that this is the remnant of the Cb seen at 22:00 GMT. The 
large increase in the area covered by cold clouds with decreasing brightness suggests that 
the storm is in its dissipating stage. Consequently, any rain received at Lethbridge should 
have fallen between passes. This is confirmed by the rain data. 

Lethbridge received 10.4mm of rain from 21-22 GMT. The weather reports, 
given in Table C.5, show that a light shower was first reported at 21:14 GMT, and 
precipitation continued until 22:44 GMT. Although no heavy rain was reported, it must of 
occured between 21:14 and 22:00 GMT, probably towards the latter half of this period, 
as heavy showers surrounded the station at 22 GMT. Unfortunately, the Lethbridge 
observer did not indicate the shower's direction of movement. 

Estimation of cloud-top temperature and brightness over Lethbridge requires 
interpolation between the two satellite passes, to correct for the storm's evolution. This 
cannot be done accurately, as the development is not a linear process, nor is it known 
exactly when the rain occured. However, some "educated guesses” can be made. 

From Figure 4.16, the temperature gradient of the cell does not flatten out until the 
-40°C contour. On the temperature map for 22:40 GMT, this occurs in the range of 


-50°C to -55°C. Assuming that the precipitation is not associated with the tight 
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temperature gradient on the boundary of the cloud, then the cloud-top temperature was 
most likely in the range -40°C to -55°C when rain was observed. As noted previously, the 
rate of change of cloud-top temperature is more rapid in the early stages of the cell's 
development. As the heavy precipitation fell after the first pass, the associated 
temperature should be a good deal lower than -40°C, because the storm was still growing 
and the temperature decreasing. Similarly, as the rain ended before the second pass, the 
temperature over Lethbridge during the rain would be slightly higher than -55°C. Also, as 
the rain occurred closer to the time of the second pass, the temperature should be toward 
the colder half of the range given above. As arough estimate, then, the temperature over 
Lethbridge while reporting precipitation is assumed to be in the range of -47°C to -52°C. 

Brightness is even more difficult to estimate. The temperature field is relatively 
straight-forward. Cloud-top temperatures become progressively lower towards the 
centre of the cloud. The same can not be said of brightness (see Figure 4.17) as values 
fluctuate throughout the cloud area, with the highest values south of Lethbridge. Most 
likely, the brightness was in the 140-150 range (as the minimum and maximum values over 
two pictures are 130 and 160). However if the area of higher brightness south of 
Lethbridge passed over the station, then the range given above Is an underestimate. 

Even if all of the above is correct, it is still unknown if Zone 1 or 2 would have 
passed over Lethbridge. Obviously this case is of limited usefulness. Be that as it may, 
the final estimates of temperature and brightness over Lethbridge during the heavy rain are 


-47°C to -52°C and 140 to 150 respectively. 
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Figure 4.18 


SOOMB ANALYSIS 
HEIGHTS AND VORTICITY 
24002 21 AUG 1979 


21 AUGUST 1979 
S/W® TROUGH PASSAGE: 
S/W RIDGE TIME: 
CONVECTIVE DAY CATEGORY: 0 
TOTAL HAIL REPORTS: 


MEDIAN HAIL TIME: N/A 


500mb and surface maps for August 21, 1979 at 2400 GMT. 
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4.2.9 July 01 1982 

The two satellite passes for this day were at 20:30 GMT and 22:10 GMT; 
respective temperature and brightness maps are given in Figures 4.19 and 4.20. 

A very large storm is seen to be moving east-northeast over southern Alberta. 
Between passes, minimum cloud-top temperatures drop from -59°C to -63°C. Brightness 
decreases dramatically from a large area with values > 160 (maximum of 170) to only a 
small area with values > 150 with a maximum of 153. 

This storm dumped more than 66mm of rain on Lethbridge in less than one hour 
before moving eastward and giving Vauxhall 26.2mm from 20-21 GMT and Brooks 
21.5mm in the next hour. As seen on Figure 4.19, the temperature and brightness over 
Vauxhall are -49°C and 157, but, less than 30km away, along the storm's track, the 
temperature and brightness are -58°C and 168, respectively. It is assumed that Zone 1 of 
the cell passed over Vauxhall and that the latter quoted temperature and brightness values 
are associated with the rain. From Figure 4.20 Brooks Is seen to be directly below the 
storm. The temperature and brightness are -58°C and 147, though brightness values of 
150 are just to the east. 

This is a straight-forward example of a large well-defined storm. Its velocity is 
approximately 55km/hr toward the east-northeast. Considering the large rainfall at 
Lethbridge and the fact that point accumulations of over 20mm were produced for at least 
2¥2 hours, this must be judged as a severe storm. 

It is also of interest to note that Calgary, on the 20:30 GNT satellite image, is under 
the storm's anvil in an area that is both cold (-56°C) and bright (162), but no rain was 
observed. This illustrates that temperature and brightness alone are insufficient 


predictors of rainfall. 
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4.2.10 July 10 1982 

The satellite passes were at 20:30 GMT and 22:10 GMT. Coronation received 
11.0mm of precipitation from 22-23 GMT. The surface weather observations, given in 
Table C.6, show that heavy rain began at 22:40 GMT and subsided by 23:00 GMT. Figure 
4.21, the temperature contour map at 22:10 GMT, shows that there are two cells west of 
Coronation, with minimum temperatures of -55°C and -41°C. As it is not possible to 
determine which of these cells gave Coronation its rain, no further consideration is given 


to this case. 
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Temperature contour map for July 10, 1982 at 2210 GMT. 
Shading indicates areas of T < -40°C 
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4.2.11 July 14 1982 

At 21:20 GMT, many Cb's covered Alberta and Saskatchewan as shown on the 
temperature and brightness maps given in Figures 4.22 and 4.23. Those in Alberta are 
embedded in a large cloud mass. The second pass, at 22:55 GMT, viewed only western 
Alberta, as can be seen from temperature and brightness maps given in Figure 4.24. 
There is little Continuity in cloud pattern between passes, in that individual cloud areas 
cannot be followed from one image to the next. However, it is known that cloud 
continued to move into northwest Alberta from British Columbia, resulting in 75-100mm 
of rain over the Smoky River Basin, north-west of Grande Prairie. 

Underneath this cloud mass, heavy rain fell at both Grande Prairie and Edson. 
Weather reports for the two stations are given in Table C.7. At Edson, continuous rain 
and stratiform clouds dominated the weather, indicating a synoptic-scale disturbance. 
Grande Prairie, on the other hand, experienced lightning and thunder with moderate or 
heavy rain reported for almost 7 hours, beginning at 15:30 GMT. Most probably, the rain 
was due to a combination of a synoptic disturbance with embedded Cb's, resulting in 
continuous rain of varying intensity. 

Since the rain at Grande Prairie and Edson is not the result of convective activity 
alone, these two cases are not suitable for inclusion in this study. From Figure 4.23, the 
cloud over Edson was neither very cold nor very bright. This indicates that the 
temperature and brightness threshold values used to infer heavy rain will be quite different 
for convective and synoptic-scale situations. 

There were 3 other locations which recorded > 10mm of precipitation on this day. 
Vegreville received 12.5mm from 20-21 GMT. Figure 4.23 shows a small bright cloud 
slightly north-northeast of Vegreville. The cloud-top temperature and brightness are 
-59°C and 164 respectively, although over Vegreville the values are only -47°C and 160. 
The temperature gradient is very tight. After shifting the cloud by 10km (the uncertainty in 
position) the cloud tops over Vegreville could be as cold as -55°C, with brightness 
essentially unchanged. The movement of the cloud is unknown. However, a general 
easterly motion would place Vegreville under Zone 1 or 2 of the cloud, with the latter 


more likely. 
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In Saskatchewan, Outlook is seen to be under Zone 1 of a large cell with a 
temperature and brightness of -60°C and 154, respectively. Brightness values over 160 
are only slightly to the west. Unfortunately, the hourly rain data are missing. The 24 hour 
rainfall total was 24.0mm. Slightly to the southwest, the same storm gave Elbow 23.8mm 
from 22-23 GMT and 24.2mm over 24 hours. Given this situation, it seems reasonable to 
assume that the 24 hour accumulation of 24.0mm recorded at Outlook was also due solely 
to the passage of this one storm. Elbow would have passed under Zone 1 of the storm 
and, if no drastic evolution took place, the temperature during rainfall would have been 
-55°C to -58°C with brightness in the range 155-170. 

Saskatoon Airport, crossed by Zone 2 of the storm, received 5.4mm of rain 
between 22-23 GMT. The temperature and brightness associated with this rain would be 
-55°C to -57°C and 155-160 respectively. Nearby Saskatoon Water received 4.8mm in 
the same time period, and with the same temperature and brightness values. 

Further south, a small cell is seen to lle over Moose Jaw. The temperature and 
brightness overhead are -54°C and 145, respectively. However, the extreme temperature 
and brightness values are only 10km to the south-west, at -58°C and 165. The cloud is 
very small and were it shifted northwards by 10km, it would miss Moose Jaw entirely. 

Only 0.7mm of rain was observed from 21-22 GMT. The weather observations, 
given in Table C.7, report light rain-showers from 20:12 GMT to 23:00 GMT with a 
thundershower at 21:17 GMT. Since there is no special report noting the end of the 
thunderstorm, it is assumed that this is an oversight on the observer's part and not simply 
a missing report. Mammatus, a sign of intense convective activity, was observed with this 
storm as well. 

There is no doubt that the cloud over Moose Jaw was a Cb. As little rain was 
reported, this is an example of a cold, bright cloud that did not produce any significant 
rain, at least not at the time it was in the satellite's field of view. While it is possible that 
heavy rain associated with the storm just missed the station, this does not seem to be the 
case as the observations give no indication of heavy showers nearby. It is noted, 
however, that only 20km to the west, Caron received 6.0mm on this day. No other 


nearby station reported more than 2.0mm, consistent with the small size of the storm. 
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4.2.12 July 16 1982 


Figure 4.25 gives the temperature and brightness maps for the 20:55 GMT 
satellite pass. A large, cold, bright cloud consisting of two cells is seen to cover much of 
southeastern Saskatchewan. This cloud moved northeast and gave 75-100mm of rain 
near Duck Mountain, in eastern Manitoba. 

TBRG data are available for Regina and Indian Head PFRA, both of which would be 
under Zone 2 of the storm. 

Regina received 2.6mm from 20-21 GMT. The weather observations, given in 
Table C.8, show that a moderate shower occurred at 20:37 GMT. At 21:00 GMT, heavy 
showers were observed north to southeast of the station. Indian Head received 6.9mm 
from 20-21 GMT and another 10.4mm from 21-22 GMT. However, no accurate timing 
of the precipitation is possible and it is not known if the 17.3mm total is the result of one 
or two periods of rain Consequently, Indian Head is not included in the sample. The 
temperature and brightness over Regina are -55°C and 162, respectively. It should be 
noted also that a northeasterly cloud motion would have taken Zone 1 of the western cell 


over Indian Head and could have produced the heavy rain from 21-22 GMT. 
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4.2.13 July 18 1982 


Temperature and brighness contour maps for the first pass, at 20:35 GMT, are 
given in Figure 4.26 and 4.27. A bright cell lies just northeast of Cluff Lake, in 
northwestern Saskatchewan. The temperature and brightness figures over Cluff Lake are 
only -29°C and 141, but 10km to the north-east the respective values are -49°C and 152. 
The centre of the cloud is 35km northeast of the station. 

Temperature and brightness maps for the 22:15 GMT pass are shown in Figure 
4.28 and 4.29 The cell near Cluff Lake has moved northeast at approximately 6O0km/hr to 
lie south of Uranium City. Its temperature and brightness have both decreased. Another 
cell now lies just to the northwest of Cluff Lake with a temperature and brightness of 
-54°C and 148, respectively. 

Cluff Lake received 10.0mm of rain between 19-20 GMT. This is consistent with 
precipitation falling from the cell shown on the 20:35 contour maps. It is also noted that 
2.6mm were received from 21-22 GMT. This might be related to the cloud north of Cluff 
Lake shown in Figure 4.29. However, this second period has not been included in the 
sample. | 

Other stations in the area reported less than 10mm. However, North Battleford 
and Cold Lake were in the path of large cells. It is seen on Figure 4.26, that Cold Lake lies 
under Zone 1 of a large cloud mass. The temperature and brightness within 10km are 
-43°C and 154, respectively. Colder and brighter values (-47°C, 158) prevail 35km to the 
northwest. This cell moved northeast at 55km/hr and, as seen in Figure 4.28, 
temperatures are only slightly higher at the time of the second pass, while the brightness 
has decreased to 120. 

The hourly rain amounts are missing for this day. The 24-hour total was 16.2mm, 
but how much of this amount came from this cell is unkown. The weather observations 
for 19-23 GMT are given in Table C.9. Outside of this period, only a few light rain 
showers were reported. It is seen that moderate or heavy rain occurred for almost the 
entire hour from 20-21 GMT. Using 5mm/hr as an average rate for the moderate rain, 
and 10mm /hr for the heavy rain leads to an estimate of 8mm during this hour. This is 
likely an underestimate as heavy rain during convective showers is often much greater than 


10mm/hr. The only other period of moderate rain was 19:33-20:00 GMT. Assuming 


73 


ee ous HHID yova a4 2 . 

Ser bes BS ate Ra eee a br} tene-nirice aft aed 
noase ant to resonthon még. ot buote ett 
Swe ri ie ae et SSO ee eqns! 
gf at lS Wlotrina ade aa vio) 63), Sail tO man Nos eT GSH - 

“MQod dpossetasn tod overt ce in bet y abdadatiet gc ‘Vid mien NO a 
lo aaadirignd G08, eUvEWwR TS? bs » bi sn ee * paewetion ont of teuj sell , > 
3 visuncaaaey Bet bre 2 
rdiwteateenod 2) 2aT TMG) OSE) “envi rier Joinwnl O} pews ete) Muld, ) 
tet baron Azle ai th 2a quoTno) aha Rar no aysors: tis set mort griliet nensiaioag 


Ps 
ae 


7 


FOND teeApannunt> ato) bafaia 7 Sa 7G" of aw... TM SS-f S ot hevisger srew mma. S 
ai AY bebo! Hens Tor eer! Doviag Snioee 2irtr savewoH -€&.5 omg a w Madi exet . 
algerie 


nap eed tral! JevswwoH anna) agrt 2241 hartge) 68% ita arouse tO 
VPA no coset! alte ents to itaq er wi bie tiled ho bret 
neern cuals sow 4¢ te Pf anoS yebad 


exis tlad.er! OS 
ene! Atiw eeeraighe bas swie egy aft) 
aff ot GE Neyer (24 °C D>] ceuley wingd' brs WSR “y'qtioedeey bet bee DEB 
BS say pl rigae ee orig TASKS Jer tear wore gevorn laa ait teswitnon 2. 
gdartidignd sairsiiiw e2acd mncone Sy 1S MT? afin -Sctpit @igiz yiete ee 

wi 0%! ot basserob earl 


onset dew lated wee BS ant yb -etelt-as? qitimdicn etme atymdonen set yiwort eT 
srolevnesdo caeissw atT nwo © beg ae ote + auneg. MWoOIne ait +o rowers wor Aad % 
(is Jig Wate yito \Woleg ev To eylatud (GO plant A mavig ae TM es-et rot 
aat roma 702 beTwIGe' mn vvese 1-8): te swe aAee al -heaQer- enw newort 
pe’ Sterebur Orlh 107 ete) spenevG Oh cu VIVIINE Orel TM 1S-OS mont wort jar 
2 ant swort alt eninge THB Ano stemitas oh Ot @jpmel ches corned a 
nacit wine rg fours! to giaraworesvuoeynos Bub nisl yvaer a6 ie 
povrmaad: TMS OOMSCE A) ‘cew niey stswoony 40 Honea como ving ait 9 


Pe 


+ 
i 


i a) 2 


2.5mm of rain in this period and 2.5mm during the rest of the day (excluding 20-21 GMT) 
implies 11mm may have fallen from 20-21 GMT. Accordingly, the rainfall total for the 
hour is estimated to be 8-1 1mm. 

One of the coldest cells seen on the 22:15 GMT image lies east of North 
Battleford with cloud-top temperatures as low as -58°C. The cell is not very bright, 
however, with only a few pixel values greater than 140. One bright patch is only 20km 
east of North Battleford and may have passed over the station. It is not certain where this 
cell originated. From the 20:35 GMT contour maps it would appear to be due to the 
intensification of the small cell that is 35km west of North Battleford. Assuming this to be 
true, then the bright patch would have passed over North Battleford at 21:20 GMT with an 
estimated temperature and brightness of -51°C +3°C and 145-150, respectively. This 
area is quite small, less than 10km in diameter, which is of the same order as positioning 
uncertainties. Consequently, either Zone 1 or 2 of the cloud could have passed over the 
rain gauge. North Battleford received 5.8mm from 21-22 GMT with 5.0mm in only 5 
minutes. From the weather reports, given in Table C.9, the heaviest precipitation 
occurred at 21:32 GMT, very close to the estimated time that the bright area passed over 


the station. 
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Figure 4.27 ~~ As in Figure 4.26 except over northern Saskatchewan. 
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4.2.14 August 11 1982 


The temperature and brightness maps for the first satellite pass at 20:50 GMT are 
given in Figure 4.30. Numerous cells are evident over central Alberta. Rocky Mountain 
House lies under Zone 1 of one of these clouds, with overhead temperature and 
brightness values of -49°C and 147, respectively. Minimum temperatures of -53°C lie 
15km to the north. 

Considerable change has occured between the two satellite passes, as seen from 
the temperature and brightness maps for 22:30 GMT, displayed in Figure 4.31. The anvils 
of some cells have merged, while others have decayed and are no nee visible. Overall, 
brightness over central Alberta has decreased. The area of temperatures < -50°C that 
was near Rocky Mountain House earlier has moved east-southeast at 45km/hr to north of, 
and over Lacombe. 

Rocky Mountain House received 9.7mm of rain from 20-21 GMT. The weather 
reports, given in Table C.10, show that the heaviest precipitation began at 20:20 GMT, 
which is 30 minutes before the satellite pass. After adjusting the cloud position for this 
time difference, the temperature and brightness over the station is -47°C and 145 
respectively. These are the final values adopted as associated with the rainfall, without 
any correction for possible evolution of the cloud. The temperature and brightness over 
Lacombe, from Figure 4.31 are -50°C and 134 respectively. An area of higher brightness 
lies south of Lacombe, but there is no evidence it passed over the station. Because of the 
great amount of cirrus debris it is not possible to discern the location of the convective 
core. Based on extrapolation from the first satellite pass, either Zone 1 or 2 of the Cb 
may have passed over Lacombe. 

The hourly rain data for Lacombe shows that 5.4mm fell from 22-23 GMT and a 
further 4.0mm from 23-24 GMT, with a total of 7.7mm in less than one hour. It would 
appear that the same cell that gave the heavy rain to Rocky Mountain House is also 
responsible, in part, for the rain at Lacombe. However, the area of temperatures less than 
-50°C has almost passed Lacombe at 22:30 GMT and yet 4.0mm was reported from 
23-24 GMT. It may be that the small, cold, bright cloud near Rocky Mountain House at 
22:30 GMT intensified while moving eastward and resutled in some of the 23-24 GMT 


precipitation at Lacombe. The amount ascribed to the 23-24 GMT period is, therefore, 
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5.4-7.7mm. 


4.3 Summary 

A summary table for all cases is given in Table 4.2. This includes cloud-top 
temperatures and brightness values, precipitation amounts and the cloud zones which 
passed over the station. Table 4.3. presents the duration and intensity of precipitation 


where this could be obtained from the TAB31 charts. 


Table 4.2 Kaine ccuml lations, 
Cloud Top Temperature (T) and Brightness (B) 
For Case Studies Discussed in Chapter 4 


Date Statwon Rain T B Zone 
(mm ) ee Co (counts) 

07/07/79 Edson tales! a) 170 1 
10/07/79 Edson ities} -49 156 1 
Penhold 26 -42 178 S 

Rocky Mtn House 0.8 36 = 165 3 

Swift Current A 3.9 -49 os 25 
Saskatoon VOR lay 1.65 1 

Out look Bai =a 5t5 160% tae 
11/07/79 Edmonton Os -48+8 165 1 
15/07/79 Out look lider -48+8 15525 1 
29/07/79 Calgary 1030 -48 168 1 

05/08/79 Cold Lake 24.0 =57 ae 1-2 
Edson 144 =o 0 156 1 
Rocky Mtn House Oy, OS) 55 1 

21/08/79 Lethbridge 10.4 SOG 2s a aotS 1-2 
O17 07/82 Vauxhall ] 2on2 -58 168 1 
Brooks Dasa) Sy) ene 1 

Calgary Oe =06 162 3 

147 077.32 Vegreville 125 S55) 164 2- 
Out look PY 6.0) 154 1 
E 1bow DES 8) 56 402255 163 1 

Saskatoon A ee} lens || Sires 2 

Saskatoon W 4.8 = 5641 158 2 

Moose Jaw Oe oo 165 | 

16/07/82 Regina 20 =oo 162 2 
18/07/82 Cluff Lake eh 0) -49 152 { 
Cold Lake Cy eee -43 154 1 

N Battleford 50 say ons 148 ie 
ime / oe Rocky Mtn House’ 9.7 S47 145 1 


Lacombe erg one 5) 0) 134 4 
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Table 4.3 Precipitation Amounts and Intensities 
From the Monthly Rain Gauge Summaries 


Station 
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Table 4.3 
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5. Results 


5.1 introduction 


An attempt will now be made to relate cloud-top temperature and brightness to 
precipitation, based on the case-studies given in the previous chapter. The interpretation 
of results presents its own difficulties and these too will be explained. Before 
proceeding, a warning with regard to the scope of this chapter is in order. The number of 
cases are few, the sample biased and many assumptions had to be made to estimate 
cloud-top temperatures and brightness. Consequently, it is unreasonable to expect the 
quality of data used here to reveal, in any great detail, the relationship between 
temperature or brightness to precipitation, should one exist. Instead, the goal is to 
develop some broad guidelines that can be used to distinguish heavy rain events. 

The trend of precipitation amounts versus temperature and brightness can be seen 
from Figure 5.1. This gives the average precipitation amount (rounded up to the nearest 
/’zazmm) for various ranges of temperature and brightness for "Zone 1” cases only and 
"precipitation > 10mm” cases, as taken from Table 4.2. The number in brackets, below 
the amount, is the number of storms in that range. Three storms have been excluded 
because of their large temperature uncertainty: Outlook on July 14 and 15, 1979 and 
Edmonton on July 11 1979. Also, Vauxhall and Brooks on July 01, 1982 and Elbow and 
Outlook on July 14 1982 have been kept as sssarete cases. It should be noted, though, 
that for the days given above, the precipitation at the two stations resulted from the same 
storm so that these are not independent cases. Figure 5.1a and b each contain 14 points 
with an overlap of 10 cases between the two. 

The trend is very clear. Precipitation amounts increase with decreasing 
temperature. For all cases for which there was > 10mm (from the TBRG hourly amounts) 
the temperature is < -45°C. Changes in brightness do not seem to be correlated with 
precipitation except that a low cut-off is apparent. Heavy precipitation cases are all 
associated with brightness > 145. 


These trends are now discussed in more detail, beginning with brightness. 
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5.2 Precipitation and Brightness 

A graph of precipitation amounts versus brightness for heavy rain cases is given in 
Figure 5.2 along with the least squares regression line. The correlation coefficient is only 
0.08. Brightness values were not corrected for sun angle. Also, its rate of change is 
unknown so that the values adopted for those cases where the time of the satellite pass 
and observation of rain differed may be in error. That little correlation was found 
between precipitation and brightness may well be a result of these two deficiencies. All 
that can be said is to repeat what was noted earlier, namely that there may be a lower limit 
of 145 that must be exceeded for heavy rain to occur. It is also possible that this limit is 
not constant, but may depend on the time of year. 

Attention is now turned to cloud-top temperature which, based on Figure 5.1, 


holds the best hope for success. 


5.3 Precipitation and Temperature 


5.3.1 Rain Amount and Temperature 

Precipitation amounts versus temperature for Zone 1 cases and heavy rain cases 
are given in Figures 5.3 and 5.4, respectively. Data for the locations on July 11,14 and 
15, given earlier, are excluded. Also, rain amounts and temperature for Vauxhall and 
Brooks (July 01, 1982) have been combined and averaged as have Elbow and Outlook (July 
14, 1982). 

Figure 5.3, for Zone 1 cases, exemplifies the difficulty in attempting to interpret 
results based on insufficient data. The regression line is given for this graph, and the 
correlation coefficent is 0.55, which is just significant at the 5% level. However, this 
value can be radically changed by the removal of only two points. 

The point representing Rocky Mountain House on August 05, 1979 at (-35°C, 
0.7mm) anchors the low end of the regression line. Without this point, the correlation 
coefficient drops to 0.37. On the other hand, Moose Jaw on July 14, 1982 represents 
the only example of a cold, bright cell with very little precipitation (-58°C, 0.7mm). 
Without this point, the correlation coefficient is 0.86, which is significant at the 1% level. 


It is interesting to note that none of the climate stations in the vicinity of Rocky Mountain 
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House reported more than 1.0mm, implying that the storm did not produce much rain. 
However, near Moose Jaw, 6.0mm was recorded at Caron. There is little evidence to 
indicate that a rain-producing cell just missed Rocky Mountain House, but this may have 
occurred at Moose Jaw. 

Turning to Figure 5.4, the precipitation versus temperature curve for heavy rain 
cases (excluding Cold Lake on July 14, 1982 as the rain amount was not obtained from the 
TBRG records), it is seen that that there are two clusters of points. One is near (-50°C, 
10mm) and the other is at (-58°C, 24mm). The correlation coefficient for this graph is 
0.81, significant at the 1% level. This does not, however, take into account the 
uncertainties in either the temperature or precipitation measurements. 

Based on Figure 5.4, then, given that > 10mm has occurred, it may be possible to 
infer the cloud-top temperature to within a few degrees. It is, of course, the reverse task 
that is of more interest, and for Figure 5.4 to prove useful in this regard, all cold, bright 
cells must produce rain. Stated another way, in order to have any hope of inferring heavy 
rain from cloud-top temperatures, it must be shown that the cloud top temperature 
distribution for convective storms producing heavy rain ts significantly different than the 
population distribution for all mature Cbs. Unfortunately, the population distribution is 
unknown and it is not a trivial task to determine It. 

Two options were considered here, neither of which proved feasible. First it was 
thought that the cloud-top temperature for all convective clouds viewed by the satellite on 
the days in the dependent sample could be used. However, there are two flaws that could 
lead to erroneous results: 
if some cells could be in the growing stage so that there might be a bias toward 

"warm" temperatures, and 
2. the atmospheric conditions which produced the storms in the dependent sample 
might be widespread, so that cells would not be independent of each other. 

An alternative is to assume that, on average, the cloud top of a mature Cb will 
reach the tropopause. If storms that produce heavy rain are colder than average, then the 
minimum cloud-top temperature should be lower than the average tropopause 


temperature. 


30 


Pee 


Went Bint te Seposeui m vide evi oF_9 ie Jari se? Dre Teee evr Yo ain inh 


ae . 
: | : 7 ; | 7 iy ih oe | on 
Ais f wert BINH iq: = . ae ae salt vite ere 


abana agin pa 
avast yaire eid? tod sea mi Ff es gen - 


Niet weer Net ovUD ot euarey riniynicinalielimanesiont 
att nett co ae 2 SDE¥ at vib 9 se bl2 geNbutoxe) ea ; ee 
QROB-) saan'al AOD .atniad Aad efopmule ah ore ornytt Jerk tarts rime 2 ARraoe 2 Aer. 
ei /rigany ale tot. trailor tego nbpaianies, anT ferent D881 te 2 vertn rit brik or 
ant worse Of pam TeveNNe’ DOR suet ‘tt lien ACL eet te racine 
ainenerienars athe gE NGO SEGRE 
hat ielieasae aq. vert) ea uci. coh HAO <TR nevi no bz 2 gayi no here® ; 
dat GOTT NSE Be299 tc ula) “Paige ware iia od ‘srutersgens? Gehbunls act xa! 


und yrinioth to SQ5¢ yn 2vsc OF ipbhio A Vey sartions Ver? chery Souler taunt 2 - 

pivhemermatpeed bunt dtl iat inured od Reap te zene aairer qou-buolo mot ms : 

art Sen TAM Ne vitae toad x) Pite, ivenrl Sng. She arniaey evtevnes vat Oly 

eS cARdieG POlldibets oo YiatarmeIn “Bed aura ie Wot recmueal noite 
fenittetst ot dea Idiver a tonal tibnw 

ee 1 7eVA SeldleseT beavers AaM.a0. veetien rer ho ahianes Baw anpitae owT 

act atilletee er yo vewely abuols svitevrics 1646) Siege Goh euers ort fartt tiquertt 

i 


; ¢ _— Sk B. : ar 
claos Ten) awall Swi 91s 'siedr avewsH. bee, adblucs stqmse tnebtaqet stm aysb arth 


, 


atures? cameo Tie of baal — 
La 
Siuworesid © ed iripim ster ter se spate griwene srl ni eo Blue ges Smee ih 


Dith 2a weceet “era” 

aiicnne rrabesgeb SM fl amie all veaeoy Apiniw anovibrica onenagarns, ot s 
Arto MaRe TO reper@Qaunl ac tort bie Wy elles terth oe. one qeeituey oe telkinn 

lliw GD anger & %o gal oud) srt agate no dine SBOE OF ovirarretls 6m 


ard rec Stam: Meat ieblng 276 re* \y veer adhe dernier 1! = 
seUbaogoN 4p’ avd oatn mart’ ewe! ect bier. wiuneretprter qe 


on 


“S@SEBD | BUOZ JO} BuN}e4adWa} SNS498A yUNOWe uoneyidinaig €°G aunbiy 


(Do) JUNIVYEYAdUdWAL dOL-Gno19 
09- ss GS- 05- S7- 07- GE - 


0c - 


Ol 


St 


02 


SZ 


NOW INWADOV NIVd 


(LU UW) 


io: ae 
A vWeonwnar 


7 
: a 
=" es 


Ye 


; 
a 


oe 


‘saseo ules Aneay 104 ainjyesaduuay snssan yuNoWe UOIeyIdIDaIg 


(Do) SHNLIVYAdWSAL dOL-Gno19 
SS- 0S - S7- Oo 


bg eunbi4 


Sus 


OE- 


Ol 


St 


02 


Sé 


NOILVINWNDOV NIVY 


(LU LU) 


7 


t ~ =a. ; 


7 
0a- oe Oe - 23- is? aé - Gé- 


)  AUTARBAMaT IOT-aQuOD 


peesordst yvaurtot e@auleye! ves) 2ustay itt bree nodiengiett P22 ao 


Tropopause temperatures are available only for Alberta at Calgary, Edmonton, 
Penhold and, for the summer of 1979, Rocky Mountain House. The test is, therefore, 
restricted only to Alberta storms which produced > 10mm of rain. Tropopause 
temperatures were kindly supplied by the Alberta Research Council for the summer of 
1979 and 1982 (the same period as for the dependent sample). For 32 days on which Cb 
activity was evident on the satellite images, the mean tropopause temperature was 
-55°C+4°C. For only those days on which heavy rain was recorded it was -53°C+4°C, 
while the mean minimum cloud-top temperature of the rain-producing storms was 
-53°C+4°C. It cannot be said that the heavy rain-producing storms are any colder than the 
“average” tropopause temperature or, by extension, any colder than the average Cb. The 
obvious difficulty, though, is that it is unkown on how many of the 32 days selected there 
were heavy rainfalls. The average tropopause temperature may also be heavily biased. 

The only alternative left to determine if useful estimates of rainfall can be deduced 
from satellite imagery is to simply test against an independent sample. This is done in 
Chapter 6, as will be discussed shortly. 

Thus far, only rain accumulations have been considered. However, accumulation 


on the ground, due to the passage of a Cb, is a combination of: 


de speed of the cloud 
23 size of the precipitating area within the cloud, and 
). the storm's rainfall rate. 


If cloud-top temperature is correlated to rain amount, it must also be related to one, or all 


of the above. Each possibility is now discussed in turn. 


5.3.2 Temperature and Speed of Cell 

For a given rain rate and rain area, the rain amount is inversely proportional to the 
speed of the cloud. Cell velocities could be measured directly from the satellite images 
for 8 of the storms investigated here. As shown in Table 5.1, the mean speed is 
50+8km/hr. The uncertainty in a single measurement is given by the position error over 
two satellite images (20km) divided by the interval between images (100min.) or 12km/hr. 
The speed of the cells is, therefore, approximately constant. However, all of the storms 


in Table 5.1 are those that produced near 10mm of rain, so that very little useful 
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(fablewo. 4 Cell Speed as Determined From 
Satellite Images 


Sequential 


Location 


Edson 

Out look 
Calgary 
Cluff Lake 
Cold Lake 

N Battleford 


Rocky Mtn House 


Average Cell Speed 


Speed 
(km/hr ) 


50 


50+8 Km/hr 


Direction 
of Motion 
towards NE 
towards SE 
towards cE 
towards NE 
towards NE 


towards —E 


towards. SE 
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information has been gained here. However, as there is no physical basis to expect rain 
amount to be related to the cloud’s speed and there is no evidence here to suggest that it 


is, itis merely assumed that no intrinsic relationship exists. 


5.3.3 Precipitation and Size of Rain Area 

Recent work by Doneaud et al. (1983) and others have shown that the total rain 
volume (RV) is highy correlated to the area and duration over which rain occurred, as 
defined by the 25dbz radar echo. The latter quantity is known as the area time integral 
(ATI). The average rain rate over the storm's lifetime is given by RV/ATI and has been 
found to be approximately constant over many convective storms. The key point here is 
that only the area covered by the radar echo, and not the magnitude of the echo, is needed 
to determine the rain amount. This is consistent with the work of Lovejoy and Austin 
(1979b) who found that temperature and brightness could be used to determine the area 
over which rain was occurring, but not the instantaneous rain rate. Rain amounts could be 
estimated, though, using an average rain rate. 

It is not possible to determine the area over which rain has occurred for the 
storms investigated here as the rain data are available only for a handful of sparsley 
distributed rain gauges. Also, the size of the heavy rain area (associated with the 
convective core) cannot be estimated directly from the satellite imagery, as it cannot easily 
be separated from the cirrus anvil in which it is embedded. However, the crucial question 
to be answered is: do larger rain areas result in greater rain amounts? With the 
appropriate assumptions this question can be addressed using the available rain data. 

The duration of precipitation is given in Table 4.3. If the speed of the storm is 
known, and if it is assumed that a circular rain area passed directly over the rain gauge (ie. 
Zone 1 cases), then the diameter of the storm's precipitation area can be calculated. This 
also requires that the advection of the cell dominates over any changes in the intrinsic cell 
rain rate. According to Zawadski (1973), this is valid for time periods less than 40 
minutes. 

Precipitation amounts versus inferred diameter of the rain area for 10 Zone 1 
cases are given in Figure 5.5. If the speed of the cell was unknown, 50 km/hr was 


assumed. The correlation coefficient is only 0.14. However, this includes the two points 
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Figure 5.5 Precipitation amount versus diameter of the storm's 
rain area for Zone | cases. 
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with the largest diameter which are associated with a rain duration of 60 minutes. This is 
longer than the 40-minute limit given by Zawadski, and without these two points there 
would be evidence for increasing rain amount with increasing rain area. When these two 
points are included, the results are, once again, inconclusive. 

Care must be taken not to attempt to extend the results given here to the radar 
work mentioned previously. Only a very small area and brief time span have been 


considered, whereas the radar estimates apply over the storm's entire lifetime. 


5.3.4 Temperature and Rain Rate 

Adler et al. (1981) related cloud top temperature to updraft velocity and hence 
also to rain rate. The coldest clouds are considered to be those that have penetrated the 
stratosphere. This requires a large amount of energy, reflected by higher updraft 
velocities. This, in turn, is related to a faster conversion of water vapour to liquid water, 
and hence to a geater rain rate. 

The maximum rain rate versus temperature graph for storms producing heavy rain 
is given in Figure 5.6. A general increase of rain rate with temperature is seen, but the 
correlation coefficient is only 0.55, which is not significant at the 5% level. However, 
from the reasoning given above, the important parameter is not simply temperature, but 
the difference between cloud top and tropopause temperature. This gives a measure of 
the storm's penetration of the stratosphere. Figure 5.7 gives the graph of rain rate 
versus temperature difference (cloud top - tropopause) for Alberta, Zone 1, heavy rain 
cases. Tropopause temperatures over the rain gauge have been estimated from the 
nearest sounding. This requires some interpolation in time and space, but as tropopause 
temperature changes are slow, these should be accurate to within 2°C. The uncertainty in 
the temperature difference is, therefore, +7°C. Figure 5.7 shows that the rain rate does 
increase with increasing temperature difference. The correlation coefficient is 0.91, 
significant at the 1% level even with this few number of points. However, this does not 
take into account the temperature uncertainty, which is more than large enough to obscure 
any possible correlation. It is also noted that the temperature difference is negative for all 
but one case, indicating that most of the storms did not reach the tropopause. It is 


suggested that atmospheric stability begins to increase at some level below the 
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tropopause. The above difficulties not-withstanding, it appears that the relationship 
between cloud-top temperature and rain accumulation is due mainly to higher rain rates at 


lower temperatures, assuming some average tropopause temperature. 


5.4 Prediction of Rain Amounts 

Based on information presented in Figures 5.2 and 5.4, Table 5.2 presents 
suggested threshold values of temperature and brightness that can be used to infer heavy 
rain accumulations from satellite data. Large ranges of temperature and precipitation 
amount are used, as the uncertainty in each variable can be quite large. The threshold 
values are valid for a “typical” Cb which is raining over an area 15km to 25km in diameter 
and moving at approximately 50km/hr. It is interesting to note that the threshold values 
’ for heavy rain over centra! Alberta found by Wieler (1981) were -50°C and 145, 
essentially identical to those given in Table 5.2. 

Before the guidelines, given in Table 5.2, are tested in the next chapter, a final note 
of caution: as the threshold values given here were obtained from measurements from the 
APT data of the TIROS satellites, they should be applied only to measurements from these 
satellites. It is quite possible that other data sources, with different resolution, would 
yield different threshold values. For example, the diameter of an overshooting Cb top can 
be < 10km, which is near the resolution limit of the APT images. The minimum 
temperature associated with this area would be underestimated by the APT data, due to 
averaging with the surrounding warmer anvil temperatures. It is unknown if such a bias is 


present here. 
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6. Satellite Rain Estimates 


6.1 Introduction 

The best test of any predictive scheme is its application to an independent data set. 
In this chapter, rain amounts will be inferred using only the the satellite-derived cloud-top 
temperature and brightness. It should be noted that Table 5.2 actually “predicts” cloud-top 
temperature and brightness, given some quantity of rain > 10mm. Here, Table 5.2 is 
being used in reverse, that is, to preform the more interesting task of inferring rainfall, 
which is not strictly valid. Consequently, some failures were both expected and 
encountered. Selection of the sample is discussed first, followed by the construction and 


interpretation of the satellite-derived rain maps. 


6.2 Selection of Sample 

HRPT images for the period May-Sept. 1982 were inspected for days exhibiting 
convective cells with cloud top temperatures < -40°C (as determined using the grey scale). 
High resolution, rather than APT images were used because of their more detailed and 
easily interpreted grey scale. Only four shades of grey are available for the APT 
products. However, temperature and brightness contour maps, on which the rain 
estimates are based, are from the APT data. 

Clouds which were very small (ie. the size of the Edmonton and Outlook storms) 
were excluded as it was thought that the temperature and brightness fields would be 
poorly resolved on the APT images. Even though the data was restricted to only one year, 
there was no shortage of suitable cases, in fact, too many to be studied here. Those 
selected represent a cross section of situations, including small isolated cells, merging 


clouds and one very large, intense storm. 
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6.3 Rain Data 


6.3.1 Climate Data 

The main source of rain data is the network of climate stations shown in Figure 
3.2. These give 24 hour rain accumulations based on a period from 08:00 to 08:00 local 
time. Hourly rain data are available only on those rare occasions when a cell passed over a 
station equipped with a TBRG. As there are no more than 2 satellite passes separated by 
100 minutes, the climate data do not allow direct comparison of satellite-inferred to 
observed rain amounts. Instead, the satellite-derived amount represents only a lower limit 
to the observed 24 hour accumulation. Even when the estimated and observed amounts 
agree, it is still possible that the observed rain occurred outside the time period of satellite 
coverage. However, should there bé a number of precipitation reports confined to the 
path of a storm (as viewed by the satellite) then it may be safely assumed that the rain !s 
due to that storm. 

Table 6.1 outlines the conditions under which a comparison of estimated to 
observed rainfall amounts is considered to be a correct, neutral, or incorrect prediction. 
Neutral, as used here, means that the prediction is not incorrect, but cannot be proven 
correct either. Of course, wherever possible, hourly rain amounts are used as a primary 


source of information in preference to climate data. 


6.3.2 Radar Data 

For a few cases, radar-derived rain accumulation were obtained from the Alberta 
Research Council radar at Penhold. Accumulations are calculated using an average rain 
rate over 5 to 7 minutes. Rain totals on the ground are inaccurate at the beginning and end 
of the accumulation period as the storm would not have moved entirely over a point on the 
ground. It is difficult, therefore, to know if a given rain amount is due to the cut-off of 
the accumulation period, or a weakening storm. To aid in resolving this problem, rain 


intensity maps for the beginning and end of the accumulation period are also included. 
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6.4 Construction of Rain Maps 
The following describes the technique used to infer point rain amounts. Satellite 


derived rain maps were constructed before obtaining observed rain data in order to avoid 


biasing the inferred amounts. 


6.4.1 Size and Placement of Rain Area 

Not only must the amount of rain be inferred, but also the size and position of the 
precipitating area within the cloud. The size cannot be determined directly from the 
satellite images. Instead, it is merely assumed to be in the range 15-25km. These are the 
values obtained from the duration of precipitation for storms in the dependent sample, 
assuming a circular rain area moved directly over the station. 

For the dependent sample, the cloud area defined as "Zone 1" was quite successful 
in isolating the portion of the cloud mass where most of the precipitation occurred. 
However, as Zone 1 constituted 20% of the cloud mass, it was usually much larger than 
the precipitating element. Consequently, while it is often not difficult to define Zone 1, 
the inferred position of the small rain area within the zone can be quite poor, especially if 
the cloud area is very large. As explained in Section 3.5.1, the convective core is 
generally near the upwind temperature gradient, and this was the main guideline used in 


determining the position of the precipitating area. 


6.4.2 Speed of Cloud 

As discussed in the previous chapter, point accumulations are dependent upon the 
cell's speed. Table 5.2 was constructed for clouds moving at approximately 5Okm/hr. 
Inferred rain amounts can easily be adjusted for different speeds. However, as the speed 
of the clouds in the independent sample used here (determined from satellite images or 
700mb maps), were found to be in the range 40-60km/hr, no adjustments were 


necessary. 
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6.4.3 Construction of Rain Amounts 

The construction of a rainfall map based on a single satellite image was a 
straight-forward procedure. For each cloud considered, temperature and brightness 
contour maps were produced. The cloud was then divided into the three zones. A 10 or 
20mm rain contour was drawn near the upwind temperature gradient in Zone 1 if 
temperature and brightness values exceeded the threshold values given in Table 5.2. The 
diameter of the rain area drawn was in the range 15-25km with smaller areas 
corresponding to smaller warmer clouds. This is a subjective procedure used to indicate 
a growing or shrinking rain area when, over two sequential images, a storm appeared to 
be intensifying or dissipating. When a 20mm contour was drawn, a 10mm contour was 
sketched around it. No attempt was made to identify the position of the Omm contour as 
this was not investigated in the dependent sample. 

A more difficult task was interpolating rain areas between two sequential images. 
It may be that rain is implied on only one, or between images. In these cases, a linear 
interpolation of the extreme values of temperature and brightness was used to determine 
the position at which to start or end the rain contours. There are four conditions that can 
be encountered: rain indicated on both images, rain on first image only, rain on second 
image only, and rain between images. 

When rain was inferred on both images, then the two precipitation maps were 
over-laid and the precipitation areas linked together. It is possible that rain actually 
occurred only at the time of the satellite passes and not in-between, although more likely 
the cloud produces rain continually, but in varying amounts during its lifetime. 

For all other conditions given above, the procedure was to first plot the path of 
the storm. Extreme values of temperature and brightness from each image were used to 
interpolate linearally the temperature and brightness values along that path. Rain contours 
were drawn wherever the threshold values given in Table 5.2 were exceeded. This 
procedure is also used to infer changes of rain totals from 10mm to 20mm and 
vice-versa. The evolution of temperature and brightness is probably not linear over 100 
minutes. However, considering the errors introduced by uncertainties in the temperature 
and brightness values, the uncertainty in the position of the convective core, and the fact 


that actual evolution rates are unknown, it would be pointless to attempt anything more 
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sophisticated than a simple linear interpolation. 

The ends of the rain contours are left open when rain is inferred to be occuring at 
the time of the satellite pass. This indicates that the beginning or end location of the rain is 
unknown. Finally, if two storms pass over the same area, the rainfall from each is 


summed together. In this manner, 30mm or 40mm contours are possible. 


6.5 Interpretation of Rain Maps 

As the error in the position of the convective core can be quite large, (>30km for 
a large cloud) emphasis is placed on the correct prediction of rain amount. That is, if the 
inferred rain amounts and contour pattern are consistent with the observed amounts and 
pattern, then this is considered a “success”, even if the position of the inferred and 
observed amounts are not coincident. There is no guaranteed sucess since, as previously 
mentioned, it is unknown when the rain given in the climate data occurred. However, 
should rain be predicted and none observed, then this is counted a definite failure. 

If the only error is the placement of the convective core, it should be possible to 
shift the inferred contours so that they match the observed rain. This procedure is to be 
avoided unless there are observed amounts nearly equal to those predicted. When 
observed rain amounts are less than predicted, it is always possible to shift the contours 
so that they lie between, rather than over climate stations. In fact, as will be seen, more 
often than not storms pass between stations. This is expected, as the inferred diameter 
of the rain area is < 25km while the distance between neighbouring climate stations can be 
much larger. Consequently, all contours are left as originally drawn. 

The following section presents satellite-derived precipitation maps for 9 days, 
with two or more storms investigated each day. Also given for each day are the 
cloud-top temperature maps. Contours are drawn only for those areas colder than -30°C, 
at 10°C intervals. The approximate scale for all temperature maps is 1 cm = 35 Km. The 
scale for the precipitation maps are given individually. For all maps (temperature and 
precipitation) an arrow is drawn pointing northwards. Each case is discussed individually 
with a subjective rating of how well the contours match the observed rain amounts. No 
objective measure is possible since the rain amounts over the time period covered by the 


satellite passes are unknown. 
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6.6 Case-Studies 


6.6.1 June 02 1982 

The cloud-top temperature contours at 21:10 GMT are given in Figure 6.1a. A 
band of cells are aligned along the mountains in northwestern Alberta. Figure 6.1b shows 
the temperature contours at 22:55 GMT. The band has moved eastward and rotated 
slightly so that it is aligned more east-west. 

Three cells, marked A, B and C on Figure 6.1 have been identified, each with 
cloud-tops colder than -50°C. Brightness values greater than 145 are found only on the 
first image, implying that the rain ended before the second satellite pass. The 
corresponding rain contours and observed amounts are shown in Figure 6.2 There are no 
observations of amounts greater than 10mm, but three stations (Nose Mountain, Kakwa 
and Smoky) enclosed by the 10mm contour associated with Cell A are very close to this 
value. Also, as the stations are aligned along the contour it is strongly suggested thali all 
of the observed rain was due to the passage of this storm. The contours identified with 
cells B and C lie between stations so no verification of the inferred amounts is possible. 

Based on the movement of the clouds and the orientation of the rain pattern, the 
area of precipitation east of the contours at Pass Creek, Eagle and Sweathouse is due to 


some cell other than those seen on the satellite photographs. 
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6.6.2 June 04 1982 

Figure 6.3 gives the temperature maps at 20:50 GMT and 22:54 GMT 
respectively. A line of cells moved east-northeast from the mountains into central 
Alberta. Cloud-top temperatures were < -60°C with brightness > 150, although it 
decreased in both magnitude and area between satellite passes. The inferred rain 
contours are shown in Figure 6.4. A 20mm contour is drawn nearly east-west through 
central Alberta. Rain amounts greater than this are observed, but 35km north of the 
predicted area and only in central Alberta. From the TBRG records, the rain at Penhold and 
Lacombe occurred after the second satellite pass, which is consistent with a northeasterly 
cloud motion. Between 23-01 GMT Penhold received 11mm in less than an hour and 
16.9mm fell at Lacombe (1 1.3mm from 00-01 GMT). 

The fit between observed and derived rain amounts would be greatly improved if 
the contours were moved northwards. The area of maximum observed rain in central 
Alberta is 35km north of the predicted location. Generally, uncertainties of the contour 
position of this magnitude are possible due to the combination of errors in gridding, 
overlaying images and positioning the small convective core within the larger anvil. 
However, in this particular case, the errors do not appear to be larger than 20km. 
Furthermore, even after adjusting the contours, there is still too much rain inferred over 
western Alberta. 

While not a complete success, the results are encouraging. Penhold and Lacombe 
both received > 10mm in less than one hour. It therefore is reasonable to assume that 
most of the rain recorded by the upstream stations also fell within a short time period. 
The rain-estimation technique was successful in that it correctly determined that heavy rain 


occurred. 
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Figure 6.4 As in Figure 6.2 except for June 04, 1982. 
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6.6.3 June 05 1982 


Two cells of interest are seen on the satellite passes at 20:35 GMT and 22:15 
GMT. These are marked A and B on Figure 6.5. The storms moved eastward from the 
Saskatchewan-Alberta border at approximately 60km/hr. The interesting feature here is 
that the clouds are aligned along the direction of motion so that point rainfall 
accumulations are a combination of rain from each storm. 

The temperature and brightness of Cell A is such that heavy rain is presumed to 
have ended before the time of the second satellite pass. The temperature of Cell A drops 
from -50°C to -57°C so that rain amounts increase from 10 to 20mm between passes. 
The resulting precipitation map is given in Figure 6.6. Moving from west to east, seen 
first is a 10mm contour associated with the rain from Cell A. Next is a 20mm contour 
which results from the overlap of two 10mm contours (one from each cell). This is 
followed by a 30 mm contours caused by the combination of 10mm amounts from the 
Cell B and 20mm from Cell A. The eastern edge of the 30mm contour marks the position 
where rain from Cell A ceased. 

As with the previous case, rain amounts are reasonable but the position of the 
contours are, at first glance, poor. Kelfield received 20mm while 39.6mm fell at Biggar, 
but the centre of 30mm contour is 40km from Biggar. Considering that only a rough 
linear interpolation is used to determine the position at which to start and end the 
precipitation areas, this is actually fairly close. Only a small clockwise rotation of the 
contour pattern is needed to greatly improve the fit between observed and inferred 
amounts. As always, it must be noted that the observed rain is not necessarily due solely 
to the cells described above. However, the proximity of observed and predicted rain 


areas strongly suggests that this is the case. 
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6.6.4 June 14 1982 

The cloud top temperature contours for the satellite passes at 20:30 GMT and 
22:10 GMT are shown in Figure 6.7. Two cells, marked A and B have been investigated. 
Cell A is southwest of Edmonton and Cell B is near Coronation. Both are moving to the 
southeast. Cloud-top temperatures are in the range -58°C to -60°C, so that the rain 
contours, shown in Figure 6.8, exhibit a 20mm contour, albeit quite small. Unfortunately, 
both cells passed over an area devoid of climate stations. The largest reported rainfalls 
are 6.0mm and 7.0mm at Alliance and Brownfield, respectively. This rain could have been 
caused by either of the cells or a combination of both. 

Cell A was just within range of the radar at Penhold. The radar accumulation map 
for 20:06 GMT to 21:46 GMT is given in Figure 6.9c. Rainfall amounts in the range 
10-25mm are common. The rain rate maps, given in Figure 6.Ga-b show that the rain 
intensity was also heavy at the above times. 

Again, the satellite-inferred rain amounts are consistent with the observations so 
that this is scored as another success, at least for Cell A. It is noted, though, that the 
radar echoes place the track of the storm east of Camrose, rather than west, as shown on 


Figure 6.8. 
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Figure 6.8 As in Figure 6.2 except for June 14, 1982. 
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6.6.5 June 15 1982 


On this day there were numerous cells in Saskatchewan as seen on the infrared 
images at 20:20 GMT and 22:00 GMT shown in Figure 6.10. The motion of the clouds 
was towards the southeast. Minimum cloud-top temperatures are in the range -60°C to 
-62°C, while brightness decreased from 165 to 150 between the satellite passes. On 
Figure 6.10a, the cells are quite small. By the time of the second pass, they have grown in 
size and the cirrus anvils have merged. Because of this merging, it is very difficult to 
follow any particular cell from one image to the next. Four cells, marked A-D are shown 
on Figure 6.10b, but the position of the convective cores within the cloud mass are very 
uncertain. 

The satellite-derived rain contours are shown in Figure 6.11. Amounts over 20mm 
are predicted, but the majority of the observations indicate much less than this. Rain is 
also seen to lie between contours, but this is misleading. At Bad Lake, between cells A 
and B, for example, 4.2mm was observed, but the TBRG records show that this occurred 
after the satellite passes. More disturbing is the observation at Dana which reported less 
than 3mm, but is enclosed by the 10mm contour associated with Cell D. However, some 
points can be made in defence. Cells A, B and C appeared to have skirted around climate 
stations so that there is no firm evidence that the contours are incorrect. Near Cell D, 
Colonsay received 19mm. There are four climate stations in Saskatoon. The minimum 
reported was 2.0mm, and the maximum, at Saskatoon Central Avenue, was 11.0mm._ If 
this observed rain is associated with Cell D, then the position of the rain contours are very 
poor, being displaced 30km to the east. However, as mentioned above, the position of 
convective cores are very uncertain and errors of this magnitude are quite possible. Note 
also that, the contours were shifted to the west, then Dana would be removed from the 
inferred heavy rain area. Overall, though, amounts appear to be greatly over-estimated so 
this case is considered to be a minor failure, at best. 

It is noted that Outlook received 5.6mm while, only 13km to the southeast, 
Glenside reported 25.0mm. Such large variations over a small scale would not be easily 
detected by the method employed here except under very favourable circumstances. For 
example, the rain-producing cell would have to be exceedingly well-defined on the satellite 


image, without a large, obscuring cirrus anvil. 
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6.6.6 June 21 1982 


Rainfall from three clouds has been mapped, based on the satellite imagery at 
20:50 and 22:30 GMT. These are marked A,B, and C on Figure 6.12. In southern 
Alberta, Cell A moved north of Cochrane as its cloud-top temperature dropped from 
-51°C to -55°C and the maximum brightness decreased from 160 to 120. The resulting 
inferred rain contours are given in Figure 6.13a. This is yet another example of a cloud 
moving between, rather than over, climate stations, so there is little support for the 
predicted amounts. The largest observed amount is 5.5mm, just upstream of the 
contours. The radar-derived accumulation map for 20:41 GMT to 21:43 GMT is given in 
Figure 6.14c. Rain intensity maps for the same times are shown in Figure 6.14a-b. The 
accumulation period was stopped mid-way between the times covered by the satellite 
passes, as the cell was clearly dying. Amounts > 10mm were detected by the radar, but 
only upstream of the satellite rain contours. Over the area covered by the contours, 
satellite estimates are greater the those of the radar. On the positive side, even though 
the cell moved over the area near Acme, no rain was predicted, because of the low 
brightness values, and none was observed. 

Two cells, marked B and C, in the vicinity of Whitecourt were also investigated. 
Precipitation was inferred to have started after the first satellite pass for both storms and, 
for Cell B, to have ended before the second pass. The resulting rain contours are given in 
Figure 6.13b 

The 10mm contour for the Cell B well represents the observed amounts. The 
22:30 GMT weather report, given in Table 6.2, shows that the storm passed directly over 
Whitecourt, yet only light rain amounts were expected as the brightness was less than 
145. This proved correct as Whitecourt received only 1.6mm from 22-23 GMT. 

The contour associated with Cell C does not match the observed amounts as well, 
but the error in position is within estimated errors. Tony Hill, which received 12.5mm, is 
less than 15km north of the contours. Ansell, south of Edson, received 14.3mm, but as 
the temperature over Ansell was higher than -10°C, the rain probably occured at some 
other time. This is supported by the Edson weather reports given in Table 6.2. During 
the time period covered by the satellite images, the only storms observed were either 


over or north of the station. Also, Edson received only 2.6mm from 21-23 GMT. 


124 


pula i tae at adaeligaenaalaae + 
ert ro} trogaun siti al ent pe aznaltete efanyia ove marth Ted .neewted grivarn, - 
art Yo risattequ, tei may st evioinenat ead raagyal sd? .emnworne betpiberq 
nb cavanah TID 28.1 S ot TMD tees yoi Ase) Mortakeruaoe tov wb-sBbe adT .2worrnoo 
eT ehr O wweit ai wore ome ranitemas, dit 10} earn yNanet nia ORT 8 owgit i 
aHilasae oi yo wersv0d eat erit megwted yaverbim heaqete 2aw bored | noltelumugos “= 
jad reper art, vd hétasteb saw mre: < eialors -oruvip vitewln eaw iso ert 2s 26988q 
sping: Sf Vd Berevon este oi! vow) Sugtags mer afiiieias ert? to meetequ ying 
roved? wave obie BViidog ant at ‘ies ant to srorll ort 14iEe ty Ov asternites stitletse : 
wal ans 49 seyeted .uetoiberg 2aw nist On creo agen aes Ort eve baviorn Rae San 
_Revwigede saw enon bag soune anentright | = 
bataoireove mole anew "wopslintW to Qpigisiv’ BAT mi 3.tre 2 team elles sw 
Bie ariete Mod 16) 2obt) Gtlisthy seul’ orig 2ST ts Ost ‘nye vA OG betwin sew oisiqued 
avnovia sxe 2Wornoo rie HomMves) a “eabG bongs art e7o'lsc babne ever or 6 led yr 
. dé! 8 ewer 
sqF 2inwore bavieeds art etreceiqay lipy. 8 Na af wo? sores ram)! @iT oe + 
Tave Vi ial haskeq moeise orlt tant aw ah SG phteT ni nsvin Sogn GL — O&;S¢ 
nar eel 2aw deeinebiid get we etsethe stew 2tweme nist trip yine May uoaetinwW 
TMO E¢-S8 mot mma.) ylac apviEsaTwosettW ts fseney beaver etl det 
llaw26 ainwone beviwede arti rose fon sack 2 Ws ritiew baMDoees wore. sdT - 
> whe. S) Beviapetirisitw 4K vooT .cias beiamitell nicziw 2 noitiaag i yon welt ud 
20 tod Une. BY beWeeey oot? Yo roe Wear. -enuerIoS ee to hoe matt mart 2aqt 7 
amuse t bex50 yidedorg, nist ert Sl eee Wann FevO: erutevegns? ant - 
eriwe aed annoy Se aver sbtpshcogs ro i 
/ [So =a 
ee oe iiioanl papa - 


125 


: Hi4urs 
“IWS O€ZZ (Gq Pue LWID OGOZ (e 1€ ZBEL ‘1 z euNr 404 yda0xe |g sunbiy ul sy 
(q) 


AseBies 


Ase6Blesy 


e 
PLouueg 


PlLouueg 


Oo . 


O€-/ (>) 


VINODOLIUM 


a 
SJ 


(3UT Uoyuowp3s 


719 aunbi4 
(e) 


. 
uosp3 


. 
PINGIAYtL UM 


isni Apineees 


o\ ganes 


— 


The) OE CS (4 bus [MD 92a0* hE 
|) OES S (4 bas TMD D208 tes -Se ad : a c 
it 0 is is COE) FS aad tol iasoxs .o oeplinié 


i fel) | . 
| 2) 


A Si#souprA 


0.0 
acl ° Bowden 
Sundre RS 
0.0 
° Olds 
505) 
°Mockingbird LO 
+ 
elas as ° Madden 
Ghost RS Ay 
Keystone LO [SS 
10 
ORO Acme CDA 909 
° Cochrane 900 
° i Carbon 
Kananaskis 
Boundary RS 
————————————— EE eT 
50 Km 


Figure 6.13a As in Figure 6.2 except for June 21, IS82: 


0.4 
°Trochu Town 


90:0 
Trochu Equity 


126 


ee Ae 


ie “4 


iy hwod : 2 ] gabe 7" 
me a 
0” ‘ 
oe 
hy wridpalatt 
. oa ’ 
dad 
5 y f 
a jot be nemtliiatt? Pai ; ic ag¢ 
; od 
7° . a —— as wot a 
pd use ‘ire 
i. OF 
1 
ALS amor 1” 
; enargoe 
Valls ~ 
Pet f abicenen 


| 2 aw 


na ee Rene ae A 


127 


“eVaqiy uiayysou 439A0 }daoxe e€ |g ainBigulsy qE|-9 esnbi4 


et roubey. 
Oat Ad Ingman 
OOM i 
SIAaUaTD | P oe : wy 0 
| z ee 
SH t2249 plod 
19 
UE ON gia 2 acs os ajepoutg 
O71 49949 Jos4e) rare Goll of! [| asuy 
01 eban, UeUOY | fh 
0°0 0°0O Ol ues par 
Aemy6 1H | BVNG 
0°0 
e6 
ac 
OT LLtH woy 01 P2eqg0 
eV Oahu hssy 74 : ° ao 

0°0° 14N0297 14uM S°ZL 07 Usaqhey, Zag 

° ° 8 * t 

. 01 Al 4eqow 
. ° 
£ qoqhey 5 4 91 


01 Atstaqatyony 
° 
Bes 


2012 


OT 490e49 Sseq 


f ae — 


07 pue| jog. 
Sho? 


01 Auoy 


eer Fi ouowIS 
O01 AIIa8UOWIS 4 


gl 


07 Axyouis 
GAG ° 


nf  “. 


nS TG, 
- (a6 a 


1 dk as ? ; 
= . ; : 
ti exe tal aie “ An Vi ,Gin suyanel ; 5.8 : 
aca od Mint vA 3) bs 
; Pa vs 
A Se 
: 7 ¥ y= 
= a _ 0.0. a PA 
yeaite ene nS Le 
8.6 0. Of Pa wices 
¢ ‘bi Bes? neon . il as 
: . Mt fer) tered Fer 94 O) bivey 7 
: joa) $07 * i455 © ade pail 4 
oa te . 
- ® 2A dauwrd big ; | : 


: Gs ie 
7 eivonsta® Pie te rain 02 a 
7 Ce . 


et ienen 7 _ 


Enadi& mailiion seve tqsaae SEG eumerni 2A des as agi 


128 


TUTELETEIIUETIUTTELEUTITET EP ET TEE EU ETE TTY 


TUTTE 


ttt 


yt 


tt 


} 


i 
tT 
yt 
tT 
t 
y 
tT 
tT 
t 
a 
T 
a 
t 
T 
} 
y 
43 
tT 
12 
T 
T 
TITT Tt 
I 
T 
Tt 
qT 
iE 
T 
a 
it 
t 
as 


T 


tT 


TTTTTITT 


TETTTETTS 


SPIO° 


Th PRETETIEREE Tt 


T 
imi 


tT 


eos 


Tt 


T 


(9) 


i! 


LIND Ev LZ (9 ‘LIND GEOZ (2 18 ZBBL ‘LZ aUNL Joy ydaoxe Gg auNbi4 UI Sy 


Freby e) 


TtTTT 


TITTY 
ttt 
TT 
tt 
Erdtyr Tt 
ELetre TY 


LESE PHY vt 


PEEEPP vt 
pritr tT 
prprr tT 

bY vt 


tT 
tT 
tT 


TETTTITTIITETT FEET TET EFT 


y 


ai § 


“LIND Evi - GEO (9 pue 


(Q) 


TEFTTTL EFT IF TEE EP ETE FET ETE ETE EL UF UF TEE ETT TTT Tt 


af 
T 
T 
af 
tT 
t 
af 
; Fielie> 
TtTtT 
TUTTI T 1 oO OO 
T TUITYT 
T TITY 
us TTT 
T fas TT 
qe SELLL rt 
T L242 yas 
T wa T 
t yore 
T hr EE 
T J 
TTY TT 
TETTETTY F FEET YITtT TT 
t T ims a 
T TTY 133 
t tT TY 
a § TT YT 
Tt eT 
T TY TT 
t SPIO* Vo apung ot 
is 1a! y 
T TT TT 
T t tT 
Tt i! vt T 
TT TEVTTET IIIT tt TJ 
Tt TT T 
tT T 1! 
A T T 
TT bn Si i Ge Si ah em) SSS Sh a GS 
tT plowed * t tT 


TY. vy = 


vi 9 eunbi4 


() 


RAPER RRO bak kok ibs ho bubs babi ap kee ib eR boo Kota ep ba bob lbon 
tT 


CiebjeD 


T 

T 

tT 

} 

t 

y 

t 

T 

TTY 
TTTTTYT af TETTTIITHST 

t THTTT 

Tt TrtTT 

af TTT 

ag tT 

tT TY 

tT VAN Ate TT 
} ELYEY TT 
t Lease 

af LLOSEe Tt 1 
tT WELEEEL T 

TT WLSSEELES Tt 

T TITETITTS 
tT 

tT 

tT 

as 

af 

tT 

tT 

y 

t 

T 


tT 
CTTTTTYT Lyvrvrdyrss TY 
TT pryryy T 
FUE FYE TT 
TIyyy Tt 
Tt Ff 
TTT 
TT TT 
T aspung °F 
TT tT 
TT TT 
tT tT 
af Tt af 
ST TIETTTICEFTT Ft T 
TT TT af 
tT t T 
t Tt T 
if 120000999 99S9R9R9 999990 0008. 


rt PIN? t 1 


SPlo: 


- 
he 


- 
nahi 
” 
* 
. 
oo 
- 
~ 
iy 
- 


CEcch rs iz 
e hb yaes ha 
tates zy 


: 
be ee be oe 
~ 
~ 


2 


ete =< 


salar 
pepcracy ; 


_ 


i= tg 


be ee ee 


7 


» a 
_ Patertcesieeettersresciirss crtsat rigrés 


ry 


fadi- 


a) : 
_ 


TMe 


Eel Sid Mp eeo 


= 
+ m 


ae ae, 7 
sagas caaeeaaes ines 


tt — 4 
ye cheptiriacd 1 _ 7 
1 r “= ‘ 
y i 
t "< : ac 
t tt ! 
) tefc- pr ate les - 
a 1% at 
“ ae 
i sia 
t er ~ at 
i ‘si + 
t i za t - 
dnteeics 2&2 etacririz: ry 
rat tl 
i t 
t sry 
£ izt 
a" r ik ‘ 
Vie str ei 
2 : ei7tt ar 
i te Pe | - 
t sat 
2 Pee ae 282s 
i ‘ ‘ j tists 
ericeizi” J taxgrrrih 
bike 1 
y i ae | perk 
t 
i 
t 
‘ 
‘ 
LaTETRCDEPELLV PPLE CA ASE CTI Lis gisbsieticogpits 


i 
Bf aie SBP’ (CC emd 107 tgao*see AO apipir 


v2ZA 
Sew tS 


oe abbpanenga : 


fame A ae ts 
~ 
- 
. 


ocberts-i *paanaueit 
mr 


= te oe be be Oe me he fF 


th: 
ij 


om Ot be te te te ee fe oe 


“pa SAGEPERUTEAEASESERIOEL UE VEee ee Ee EPETAsEs 
= 7 - 


te! 


BT 4 snupi4 
ris 


Table 


aT 
Men 
Gree) 
eek 
NET; 
eal 


6. 


Se 
RS 
Sie 


SA 


SP 


SA 
SP 
SP 
Sie 
SE 
Su 


2100 
2200 
La 


2500 


2eA0 


2100 
209 
Edge) 
2100 
Zale 
2245 


20SCT 


SUSET 


SOS]. 


A0SCT 


40SCT 


A0SCT 
A0SCT 


SURFACE WEATHER OBSERVATIONS FOR JUNE 21, 1982 


WHITECOURT 


655Cl 2o0SBKN 157 CUZCCESCIi 
So oC is eIUBKNG To sim CUCCE SACOG I GCGrow 
E60BKN 90BKN 280BKN 15TRW- CU2CB5AC1CI 1 
TOoTM @VRHD ETG ALODS 
EGOBKN SOBKN 15+ TRW=sSC2CBOACi 
TSTM MOVG E 


EGUBKN SOBAN Wl52i SCs) CRS AC SHWRS E 


EDSON 


E60BKN 110BKN 40 CF1TCUBAC1 SHWR NW 
E60VKN 40T CF1 TCU6AC1 SHWR NW 


EGUBKNE ZORRW- (CUS E1GCG wy 


E60BKN 


AORRW SA CB Crna) lean TZ 


ES0BKN S07 =CBS [STM MGVGE 


ES0BKN SOs CBS 


129 


Oo DY 
Bote |i 


fideaosva 


+8 aMUL O03 2MOUTAVAS28 AInTAAW SOARRUE 


- 


, 
es 
— —_ 
_ 7 

aw i. 


) 


; 


i 
- a7 


TAUGIIT IHW > 


sor WA OBSITI2Z2a TIe0e WOES Ae 


yoaoe OOS en UST © 


W2 ODeTI 1 (9SoagEosua T+Al YAedess Foece 
Prgrgaedasua. wale ‘Naoao wade vieOGa TIe0e OSSS Fe OSI 
eGOIA Ta OAnVO Wier 
a 2 y - : 
(Qhasoede -wit +2. piggeswiaosa ToROs Get Ae URY. 
+ OvVOM MfeT - Hf . 
. | 7 
5. onWwHe 1Od SE3 202 Tre@t UNabe Wikoes T3e0t OPES a2 USY: : 


i 


WH AWhe 1 OARUDT: 


Wd Rae. 122. 
W 

BNI¢e; AG 

4 


i> Ob. VABOTT 


C.J B80. -AsWaITOD MngOsa ORFs 92 a 


‘Ave Os. 


wNadds TI20b OOFS AE THY, 
HD i Waves iD dage gure 42 Tay 


View, 
| =i 6 
hyaT') GQaT -WATOS MaeOBa Sets oe Tey 


avoM: What eed TOE ceiaaal irs 927 


682, 08 , sk jal a + os 1 : » 
e 7 7 pe 


6.6.7 June 22 1982 

Figure 6.15a-b gives the temperature contours at 20:35 GMT and 22:15 GMT 
respectively. The storms of interest are those in eastern Alberta that moved eastward 
over Saskatchewan, marked A,B,C, and the cells marked D and E in south-western 
Saskatchewan. Just north of Cell A is a large precipitating cloud band. 

From the satellite images, it is clear that considerable changes have occurred 
between passes. Some cells have merged and new storms have formed. Similar to the 
situation outlined for June 15, the position and track of the convective cores is very 
uncertain. In fact, it appears that not all of the cells have moved in the same direction as 
reflected by the inferred rain contours given in Figure 6.16a and 6.16b. 

This is by far the most complicated of all cases investigated. Generally, results are 
acceptable in central Saskatchewan, but poor in the south. Coronation and Sullivan Lake 
both received near 8mm over 24 hours, and are very close to the 10mm contour 
associated with Cell A. However, from the TBRG records, ony 5.5mm fell at Coronation 
during the period covered by the satellite images. The weather observations for 
Coronation are given in Table 6.3. Moderate rain occurred at 22:00 GMT, associated with 
gusty winds of over 30 knots. It is quite likely that the gauge-rain amounts are 
under-estimated because of these strong winds. The narrow 20mm contour associated 
with Cell B lies near three stations, two of which (Scotfield and Sibbald) received nearly 
15mm. Further south, conditions are worse. Acadia Valley lies within the 10mm contour 
of Cell C, but did not receive any rain whatever. The nadir of all test results is found over 
southern Saskatchewan. There are a number of stations enclosed by a 10mm or 20mm 
contour for which only a trace of precipitation was recorded. There is no obvious 
explanation for the complete failure of the satellite rain estimates in this area. 

Overall, the satellite rain estimates do not provide a good fit to the observations. 
However, this is by far the most complex case. It is possible that cloud areas identified as 
being the same cell over two images are, in fact, different cells. If true, then the area 
covered by the rain contours would be much smaller. A shorter interval between satellite 
images than the 100 minutes provided by TIROS is required to follow individual cells when 
they are merging and evolving at arapid rate. It is not surprising that this example proved 


too complicated for the simple technique employed here. 
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Figure 6.15a Asin Figure 6.1 except for June 22, 1982 at 2035 GMT. 
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Figure 6.15b Asin Figure 6.15a except at 2215 GMT. 
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6.6.8 June 27 1982 


The temperature contours for the single infrared image (at 21:15 GMT) available 
for this case is given in Figure 6.17. The attempt to map rainfall from this storm was 
approached with some hesitation. It is much larger than the other storms studied and may 
represent a different scale of development. The minimum cloud-top temperature is 
-71°C, far colder than any other cell investigated. Clearly, this storm is not the typical Cb 
that is Common over Alberta and Saskatchewan. However, it was thought to be 
interesting to apply the rain-estimation technique to a variety of situations, so despite any 
possible misgivings, this case was included in the study. 

Using Figure 5.4 as guidance, indicates maximum amounts of 30mm, 
corresponding to the very low temperatures. The size of the heavy rain area was drawn 
somewhat larger than usual, as a subjective response to the size and apparent intensity of 
the storm. 

The rain contours are given in Figure 6.18 The arrow points toward the assumed 
direction of cloud motion based on the 700mb flow. 

Estimated rainfalls do not match the observed amounts very well, but these are 
quite variable. Val-Marie, within the 20mm contour, received oniy 11.0mm. However, it 
is near the western edge of the cell and positioning errors could account for the 
discrepancy. A slight eastward shift of the rain pattern would put Val-Marie between the 
10 and 20mm contours. Conditions further north cannot be explained so easily. Kincaid 
received only 2.2mm while slightly to the east 25.0mm fell on Lafleche. Generally, 
though, amounts in this area are near 10mm, consistent with the position of the 10mm 
contour. Accumulations > 20mm are common east of, and along the assumed storm 
track. However, these occurred well after the time of the satellite pass and considerable 
evolution of the storm could have taken place. Ormiston, for example, received its 26mm 
three hours after the time of the satellite image. 

While this case is not a complete failure, as heavy rain amounts were reported, it is 
far from a success. As usual, more rain data would be helpful, especially in extreme 
southern Saskatchewan. The best that can be said is that the threat of heavy rain was 
properly predicted. However, with a storm of this intensity, even a novice meteorologist 


would have foreseen the possibility of associated severe weather. 
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Figure 6.17 AsinFigure 6.1 except for June 27, 1982 at 2115 GMT. 
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6.6.9 July 29 1982 


From the temperature contour maps at 20:00 GMT and 21:40 GMT, given in Figure 
6.19, it is seen that there are many small cells that intensified between passes. This day 
was chosen for study as the majority of the cells appeared to be only borderline cases 
from the point of view of producing heavy rain amounts. As can be seen on the 
precipitation map given in Figure 6.20 the rain areas are quite small. Unfortunately, once 
again they lie mainly between, rather than over the climate stations. The largest area of 
inferred precipitation is just northwest of Edson, where a 20mm contour was drawn. 
Edson received only 4.1mm from 22-23 GMT. However, the weather reports, shown in 
Table 6.4, reveal that the storm was an intense one, giving heavy rain and hail. Also, as 
only 6/ 10 of the sky was covered by the Cb it is reasonable to assume that the storm did 
not pass directly overhead, and that amounts > 4mm could have fallen nearby. The area 
of greatest observed precipitation is near Rocky Mountain House. TBRG records show 
that this rain did not occur until 01-02 GMT and is possibly associated with the small 
10mm contour to the southwest. This area has been denoted "forming’ to indicate that a 
new cell appeared at this location on the second satellite image. New, in this case, means 
that the cell was not visible on the first image. 

The 10mm contour east of Edmonton is associated with a small storm that 
intensified as it moved eastward through the Edmonton area. Rain amounts are light. 
Woodbend received 3.0mm while nearby stations reported less than 1mm. The cell was 
also seen by the ARC radar, but, unfortunately, not until 21:34 GMT. Figure 6.2 1a the rain 
intensity map for that time, shows that rates > 10mm/hr were occuring. However, the 
storm dissipated by 22:08 GMT and the rain accumulation map for 21:34 GMT to 22:08 
GMT gives amounts less than 2.5mm. The weather reports for Edmonton International 
Airport are given in Table 6.4 and show that a heavy shower was seen east of the airport 
at 21:21 GMT. 

Overall, no verification of the inferred rain amounts is possible, as the storms 
moved between climate stations. However, the lack of even one report of accumulations 


as large as predicted suggests that the inferred rain amounts are too large. 
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6.4 SURFACE WEATHER OBSERVATIONS FOR JULY 29, 1982 
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6.7 Discussion 


Based on the cases studied some general comments can be made. 

Inferred rain amounts and the position of the contours best represented the 
observed rainfall for isolated cells. Results were expectedly disappointing for cases 
involving merging cells, as large areas of obscuring cirrus debris from the merged anvils 
made it difficut to follow any given cloud area from one image to the next. Increased 
frequency of satellite images would aid in the solution to this problem. Even for isolated 
cells, the deduced track of the cell often differed from the actual path (based on the rain 
amounts) due to a combination of three errors: 

1. gridding uncertainties (10km) 

2. errors introduced when two images were overlaid (up to 30km) 

3. uncertainty in the position of the convective core within the cirrus anvil (error 
dependent on size of cloud area) 

While precipitation amounts appear to be generally overestimated, this is to be 
expected. The rain contours drawn represent the maximum amounts associated with a 
cell passing directly overhead. As the cells studied here are generally smaller than the 
distance between climate stations, it is quite likely that, in many cases, the maximum 
rainfall occurred at some location between the stations. Radar-derived accumulations 
were obtained for three storms. Satellite estimates agreed on one occasion but were 
greater than the radar amounts for the other two. Results would have to be described as 
"mixed". 

The observed rain amounts occasionally exhioited large variations over short 
distances (see discussion for June 14 and 15 for example). It is unknown if this was 
caused by a small rain area passing directly over one station while just grazing another, or 
the result of small-scale variations within a large precipitation element. The rain estimation 
method used here is not capable of discerning the latter situation. 

One extremely intense storm was investigated. In general, intense storms are not 
expected to be well handled. The largest rain amount that can be predicted, for a storm 
moving at about 50km/hr, is 30mm (corresponding to a cloud-top temperature of -70°C). 
One-hour accumulations much greater than this are possible. For example, the storm of 


July 01, 1982 (see Chapter 4) gave 66mm to Lethbridge. Cloud-top temperatures at that 
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time are unknown, but were only -60°C one hour later. It may be that some feature other 
than temperature alone can be used to infer higher rain amounts. An example would be a 
storm exhibiting a "V" shaped cloud-top temperature structure. Prediction of the rain 
rates associated with such features will not be possible until they are first included in the 
development sample. 

Finally, the question of "false alarm ratio” (the fraction of storms that are cold and 
bright but do not give heavy rain) has not been answered by the few case-studes 
investigated here. The cells over southern Saskatchewan on June 22, 1982 and near 
Moose Jaw on July 14, 1982 (see Chapter 4) were all bright and cold but little rain was 
reported. In the case of Moose Jaw, a thunderstorm is known to have occured and it is 
possible that the storm only grazed the city. However, the balance of cases strongly 
suggest that at least the threat of heavy rainfall somewhere along the storms path can be 
succesfully identified. With better time resolution of satellite images, the crude 
interpolation of temperature and brightness required to determine the beginning and end 
of rain areas would not be required. This should lead to much more accurate estimates of 
rain areas. 

Overall, then, it is felt that "reasonable" estimates of rain amounts can be made 
using the TIROS satellite images. Complete accuracy is not obtainable as even the ground 
observed rain amounts can be in considerable error. Here, rain amounts were predicted 
only in very broad ranges. No doubt, more accurate and higher resolution satellite and rain 
data would result in a better rain estimation scheme, as would the inclusion of radar data. 


This, and other important improvements are the subject of the next and final chapter. 
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7. Conclusions and Future Work 


7.1 Conclusions 
Based on the development sample the following conclusions have been reached: 
ie Given point rainfall amounts of > 10mm, produced by the passage of a 


convective storm, the cloud top temperature is < -49°C and the brightness is 


>145. 
2. The amount of precipitation is inversely proportional to temperature. 
oF Rain intensity is proportional to the temperature difference between the cloud 


top and the tropopause. 
Result (3) is tentative as the sample size upon which it is based is small. However, it is 
potentially the most important, as the point rain accumulation, based on the storms studied 
here, rests mainly on the storm's rain rate. 

When (1) and (2) were applied to a set of satellite images (in order to infer rain 
amounts) the estimates were comparable to observed accumulations for isolated cells. 
Results were poor for more complex situations such as merging cells, because individual 
cloud areas could not be tracked over two succesive images. Inferred rain ranges were 
large and the majority of observed rain accumulations were for 24 hours, so that 
comparisons are not completly reliable. Nonetheless, results are encouraging enough to 
warrant further investigation. Given in the next section are suggestions for extensions to 
the work done here that would, it is hoped, lead to a useful satellite rain estimation 


technique for the Canadian Prairies. 


7.2 Future Work 


7.2.1 Satellite 

A basic step forward would be to use the HRPT visible and infrared images. Digital 
data for these are available from the Alberta Weather Centre in Edmonton. However, 
TBRG records may not be available until many months after the period they cover which 
requires that many months of satellite data be saved on magnetic tape until the days 


suitable for study are known. This could be a problem as tapes are not regularly archived 
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at the Weather Centre. The alternative is to produce temperature and brightness maps on 
a daily basis, however this too may not be feasible within the constraints of an operational 
weather office. Undoubtedly some compromise could be reached, such as saving only 
tapes of days which exhibit considerable convective activity. 

Increased time resolution is also essential. GOES images are available every 30 
minutes, but the resolution of the infrared images is only 16km over central Alberta. This 
is still useful as areal cloud growth could be monitored. Also relationships between 
temperature fields from co-incident GOES and TIROS images could be used to revise any 
incorrect values caused by the poor GOES resolution. At the very least, temperature 
changes could be obtained from the ’2 hourly GOES pictures. 

Duration of satellite coverage is important as well. The Weather Centre receives 
data from both the NOAA-7 and 8 satellites. At 2 passes/satellite this allows coverage 
for 5 hours. Of course, the use of GOES images makes continuous cloud monitoring 
possible. It would be interesting to determine if the correlation between satellite 
estimated and observed rain amounts increased with increasing duration of coverage. 

Brightness values should be normalized for overhead illumination. This would 
allow better comparisons between different geographical areas, time of year, and other 


techniques developed by other investigators. 


7.2.2 Radar 

The most important improvement, next to better satellite data, is the introduction 
of radar data. The high temporal resolution would ensure that a radar map would be 
available at a time very close to that of the satellite pass. It would, therefore, no longer be 
necessary to correct and estimate cloud-top temperatures and brightness for time 
differences between the satellite picture and observation of rain. The wide areal 
coverage would allow correlation between echo strength or area to temperature and 
brightness throughout the cloud area. This is especially valuable as it would help define 
the position and size of the storm's convective core in relation to the overall cloud area 
and temperature structure. Furthermore, changes in the intensity of the rain rate could be 


easily monitored. 
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7.2.3 Rain 


Even with radar data, there is still a pressing need for improved ground-truthing. It 
would be valuable to implement a dense network of high resolution rain gauges in an area 
frequented by convective storms, such as central Alberta. As a minimum, the actual 
traces from the existing TBRG network should be made available, rather than only the 
TAB31 charts which give only hourly amounts. This would allow accurate timing of the 
precipitation event, both in amount and intensity. These are available by special request 


from AES headquaters in Toronto, but only 8 months to a year after the date of the rain 


event. 


7.2.4 Statistics 

A much larger sample than was used here is required for proper statistical analysis. 
Also, all convective storms should be used in the development sample, independent of 
temperature, brightness or rainfall. Of course, this will result in a very large sample, but 
more importantly, it would avoid biasing the sample. This would eliminate one of the 
problems encountered here, namely that it is not Known if the distribution of cloud-top 
temperatures for the "heavy rain” cases is different from that for all storms. Also the use 
of all cells would lead to a scheme that could predict smaller rain amounts (O- 10mm), and 


allow an estimate to be made of the false alarm ratio for inferring heavy rain. 


7.2.5 Further Suggestions 

The shape of the cloud area or temperature structure could be used to infer the 
intensity of a storm. For example, a wedge-shaped anvil may indicate strong divergence 
aloft and hence amore intense storm (Scofield and Oliver, 1981). 

Climatology should also be used to adjust rain estimates. If a geographical area is 
known to receive enhanced or diminished rainfalls, due to some topographical feature, this 
should be reflected in the predicted rain amounts. 

Finally, all tests mentioned above or performed in this study could be applied to 
synoptic-scale rain events. A preliminary test was made for such a case that occurred on 
July 02, 1982. Considerable rainfall was found beneath an area with cloud-top 


temperatures near -20°C while almost no rain was reported under an area with extensive 
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cloud cover and temperatures near -45°C. Similar behavior was seen over Edson on July 
14, 1982 (see Chapter 3). Clearly, the procedure used to infer synoptic rain amounts may 


be quite different than those for convective storms. 


If all that is required is some means to identify areas were heavy rain is likely, then 
the simple technique presented in this study is a valuable tool. However, much work 
remains to be done. Given the experience of other investigators, there is little doubt that 
a method of estimating rain amounts, large or small, from satellites can be developed. It is 
hoped that the work presented here proves sufficiently encouraging to stimulate further 


reasearch, and that the suggestions given above will soon be acted upon. 
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Appendix A 


Infrared Calibration Equations 


The step-by-step procedure used to convert the infrared radiance counts to equivalent 
black-body temperature has been described in section 2.1.2. Here, the coefficients and 


equations needed to apply those procedures are given. 


A.1 Conversion of PRT Counts to Temperature 

The equations used to convert the satellite counts from the four platinum resistor 
thermistors into temperature in °K are given in Table A.1. TMIRP are the digital counts. 
The factors of 4 and 16 are used to convert the received 8 bit number into a 10 bit 


number. 


A.2 Conversion from Temperature to Black-body Radiation 


The black-body radiation N(T) associated with temperature T is given by 


y. 


2 
N(T) =§ 6(v,T)8(v) dv (A.1) 
~ 
where 8(y,1) is the Planck function, y is the wavenumber (cm7! ) and 6(v) is given 
by 
¢(v) 


ev) = mar (A.2) 


where V} and V2 are the wavenumbers at the limit of the response function. 


Approximations of equations A.1 and A.2 are 


N(T) = 3 (v, 7) 4(v,)dv (A.3) 
and 
oy) = = fh 
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The form of the Planck function used is 


B(v. 7) = he /(exp(C,v;/T)-1) (A.5) 


where the universal constants C1 and C2 are 
C1 = 1.1910659 x 10°° mw/(m’ sterad cm“ ) 
C2= 1.438833 cms 


Valuesof vo. Av. n and $(v.) are given in Table A.2. 


A.3 Correction for Non-Linearity of Space 
As the output of the radiometer is not perfectly linear, a correction Is applied at 
step 5 in section 2.1.2 by using a non-zero value for the radiance of space. The values 


used are -1.151 and-1.176 for TIROS-N and NOAA-7, respectively. 
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Appendix B 


Calculation of Grid Positions 


The position of townsites in terms of scan line and pixel number was determined by the 
point-by-point gridding procedure described by Reinelt et al. (1975, pp l3s5-14 1): 
Therefore, only a brief description is given here, following the summary given by Greaves 
(1981). Reference should be made to Reinelt et al. for further details. 

The geometry of the satellite's orbital path, assuming a circular orbit and spherical 
earth, is given in Figure B.1. ‘The satellite is assumed to start at the equator at $0>%0 
where 9 Is the latitude and @_ is the longtitude The satellite moves along the arc "s’ 
and scans a point on the ground at position $),9) , an arc length Y radians from the 
satellite subpoint. | 

During the time that the satellite has moved over the arc s, the earth has rotated a 


distance given by 


(s/2m) (27L/360) = ks (B. 1) 


where k=L/360 and L is the satellite's regression rate in degrees per orbit. The total 
longitudinal displacement that occurs between the time that the satellite crossed the 


equator and scanned position $),9, is, therefore, given by 


Vo9ne1-Ks (Be2) 


Spherical trigonometry is then used to determine a relationship between the arc length s 


and the position $),6,.From Reinelt et al., this is given by 


F(s) = sin(s)cos(69-6,-ks)cos($,) - cos(s)sin(i)sin($,) 


Be cos (s)sin(@9-8,-ks)cos (i )cos ($}) =0 ae 


This equation is solved by first taking its derivitive 


F'(s) = cos(s)cos(@9-6,-ks)cos($,) + ksin(s)sin(89-8,-ks)cos (61) 
Sins) Simi) simi.) “Sim s)Simien-e7 ks )cos( 1cos oq) (B.4) 
= keos(5)\cOs(Cig- Oaks) cost I cOsiany) 
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Figure B. 1 Geometry of a polar-orbiting satellite on a 
non-rotating earth. See the text for explanation 
of symbols. (After Reinelt et al., 1975). 


then, using a first guess for s, performing the iterative procedure 
Si) easels PiscEU (sa) (B.5) 


which converges to the root of equation B.3. ‘The initial guess, s,, is determined by 


re-writing equation B.3 as 


sin sgo2 _ sin(i)sin($)) - sin(69-8)-ks cos (i)cos($)) 
tan So Lr eas (0,0) -ks9)cos $7) - (B.6) 


and settings, = ¢) onthe right hand side of equation B.6. 
The number of the scan line which passes through $],8) (counting from the 


equator) is easily determine from 


# of scan lines= (s)(P)(RPM)/2z (B.7) 


where P is the orbital period in minutes and RPM is the rotation rate of the scanning 
radiometer. 
The position along the scan line from the subpoint to $,,6, is found by relatiing 


the arc length Y along the ground, given by 
siny = sin(69-8,-ks)sin(i i)cos($,) + cos(i)sin($)) (B.8) 
to the rotation of the scanning radiometer n (see Figure B.2) found from 
tann = Rsiny/(H+R(1-cosy) ) (B.9) 


where R is the radius of the earth, and H is the height of the satellite above the ground. 


The arc length n_ is then converted to a pixel along the scan line using 


(n) (DIGF) (60)/ ( (27) (RPM) ) (B.10) 


where DIGF is the digitizing frequency in values per second. 
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Satellite 


Figure B.2 Geometry of the scanning radiometer as it views the earth. 
See the text for explanation of symbols. 
(After Reinelt et al., 1975) 
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Appendix C 


Surface Weather Observations 


Following are the surface weather observations for the dependent sample discussed in 
chapter 4. The format is as follows: 
ve three letter station identifier 

type of report 

time of report (GMT) 

cloud heights (hundreths of feet) and amounts (tenths) 


visibility (statute miles) 


winds (in tens degrees and knots, included only if > 15 Kts) 


2 
3 
4 
a 
So weather and obstructions to vision 
> 
8 cloud type and opacities 

2 


remarks 
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SURFACE WEATHER OBSERVATIONS FOR JULY 07, 


1900 
120 
2000 
2005 
2040 
vanele) 
2200 


E200VC 
ESOOVC 
ESOOVC 
E200VC 
E25BKN 
E25BKN 
E30BKN 


EDSON 


6TRW=-A CB10 L1GE€G ALQDS 


15TRW- 
15TRW- 
2TRW+F 
LS TERM = 
20TRW- 
SORW- 


CB1i0 LTGCG ALQDS 
CB10 LTGCG ALQDS 
CBr sEiGeG ALQDs 
CBe 

GBS IGG 

CBO METGCG, NE 


pSies 


163 


aver , To Yul AOA 2uGLTAV 


eO0dA DIMT1 OF: 
eqp1A oc 
20GgtaA O28 


AID J A DOE 


Hs) 


O2I9T 108d 


QW hoot Bad. 


k- 


yoeos 


nner 
ocer 
anos 


2208 


“WaT QvOOs3 
" VOOR’ 
svO0E3 
g° VOOR? 
WNAESS 
wer sa pct gots 

WHOS wiatee. ones aa 


-WhTet 


CAGE. 


29 
a2 * 
Aa ° 
92 
2 
me) 


164 


Table C.2 SURFACE WEATHER OBSERVATIONS FOR YULY 10, 1979 
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SURFACE WEATHER OBSERVATIONS FOR JULY 14, 1979 
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SURFACE WEATHER OBSERVATIONS FOR AUGUST 11, 1982 
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